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measurements.  For  comparison,  the  pressure  measurements  were  made  with  the  canard  on 
and  with  the  canard  off.  From  the  pressure  measurements,  flowfield  velocity-vector, 
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I.   INTRODUCTION 

A.  SUPERMANEUVERABILITY 

In  modern  aerial  combat,  the  ability  of  an  aircraft  to  maneuver  into  the  post-stall  regime 
for  short  periods  of  time  is  crucial  to  its  survivability.  Current  fighter  aircraft,  such  as  the 
Israel  Aircraft  Industries  Lavi,  SAAB  Gripen  or  European  Fighter  Aircraft  (EFA),  employ 
a  close-coupled  canard  to  allow  for  continued  maneuvering  where  conventional  aircraft  may 
have  departed  from  controlled  flight.  This  increase  in  maneuverability  results  from  the 
favorable  interaction  of  vortices  over  a  delta  wing.  Double-delta  wings  or  leading-edge 
strakes,  such  as  those  on  the  F/A-18  or  F-16,  have  been  used  to  enhance  the  lift  in  the  same 
way.  However,  until  recently,  only  the  SAAB  Viggen  had  been  successful  at  using  a 
canard  to  maintain  lift  at  high  angles  of  attack. 

B.  CLOSE-COUPLED  CANARD 

The  advantages  of  a  close-coupled  canard  have  been  known  since  the  1960's.  It  was 
found  by  Behrbohm  [Ref.  11  that  the  combination  of  a  close-coupled  canard  and  delta- 
wing,  of  small  aspect  ratios,  has  significant  advantages  over  a  conventional  delta-wing  or 
wing/horizontal-tail  configured  aircraft.  Both  CLmax  and  the  angle  of  attack  for  CLmax  are 
increased  by  the  addition  of  a  delta-canard  to  a  delta-wing.  During  the  1970's,  an 
experiment  was  performed  by  Lacey  [Ref.  2]  to  determine  the  correct  canard  geometry  and 
location  for  maximum  lift  enhancement.  As  a  result  of  Lacey's  work,  Behrbohm's 
conclusions  concerning  the  use  of  a  delta-canard  with  a  delta- wing  were  confirmed.  It  was 
found  that  locating  such  a  canard  above,  rather  than  coplanar  with,  the  main  wing  produced 
the  most  favorable  vortex  interaction.  In  the  1980's  canard  research  continued.  Work  by 
Er-El  [Ref.  31,  Stoll  and  Koenig  [Ref.  4]  and  Calarese  [Ref.  5]  provided  insight  into 


canard/wing  vortex  interaction.  However,  most  of  their  work  was  done  at  low  to  moderate 

angles  of  attack  with  little  or  no  canard  deflection. 

At  the  Naval  Postgraduate  School,  a  series  of  experiments  has  been  conducted  to 

compare  a  close-coupled  canard  model,  designed  in  accordance  with  Lacey's  work  [Ref. 

2],  to  a  wing-alone  configuration.  The  first  tests,  conducted  by  Kersh  [Ref  6],  were  to 

o 

determine  the  forces  on  the  model  up  to  50  angle  of  attack.  During  these  tests,  the  canard 

o 

was  deflected  between  plus  and  minus  25    to  determine  which  canard  incidence  angle 

would  produced  the  maximum  lift  at  a  given  model  angle  of  attack.  It  was  determined  that 
the  maximum  lift  enhancement  of  the  canard/wing  configuration  over  the  wing-alone 

O  0 

configuration  occurred  at  22    model  angle  of  attack  with  the  canard  set  at  a  positive  7 
incidence.  At  this  angle  of  attack,  the  first  stall  (loss  of  lift)  occurred  on  the  wing  for  the 
wing-alone  case;  the  canard  vortex  seemed  to  provide  a  reattaching  mechanism.  Reference 
6  is  a  complete  discussion  of  Kersh's  work. 

C.  THESIS  OBJECTIVE 

This  thesis  was  the  second  in  the  series  of  tests  conducted  for  canard/wing  lift 
enhancement.  The  model  was  set  for  maximum  lift  enhancement  conditions  as  described 
above.  Wake  surveys  were  then  conducted  using  a  nulling  five-hole  pressure  probe.  The 
objective  was  to  gain  quantitative  data  concerning  the  total  pressure  and  velocity  profiles  at 
three  crossplane  locations,  while  the  model  was  operating  at  22°  angle  of  attack,  with  and 

without  the  canard.  By  this  approach  the  effect  of  the  canard  leading-edge  vortex  on  the 
main  wing  leading-edge  vortex  would  be  further  investigated.  From  the  pressure  data, 
velocity-vector  plots  and  total-pressure  contours  were  generated  for  a  comparison  between 
the  wing-alone  and  canard/wing  configurations.  Such  a  comparison  served  to  help  reveal 
the  enhancement  mechanism  at  an  angle  of  attack  beyond  those  investigated  previously. 


II.    BACKGROUND 

A.  VORTEX  GENERATION 

The  dominant  characteristic  of  flow  over  a  highly-swept  delta  wing  is  the  generation  of 
a  strong  leading-edge  vortex  as  shown  in  Figure  1.  These  vortices  are  the  result  of 
separated  flow  at  medium  to  high  angles  of  attack.  As  opposed  to  the  chaotic,  separated 
flow  associated  with  stall,  these  vortices  are  stable,  coherent  sources  of  high  energy  and 
generate  areas  of  low  static  pressure  resulting  in  localized  regions  of  enhanced  lift  at  the 
leading  edge. 


Primary   vortex 


Secondary   vortex 

Figure  1.  Flowfield  Over  the  Top  of  a  Delta  Wing 

The  strongest  vortices  result  from  planforms  with  sharp,  highly-swept  (>50°)  leading 

edges.  The  sharp  leading  edge  promotes  the  leading-edge  separation  necessary  for  vortex 
generation.  However,  the  combination  of  large  sweep  and  a  sharp  leading-edge  produces  a 
large  planform  with  a  low  (L/D)max  and  a  shallow  lift-curve  slope.  The  net  result  is  poor 
range  and  endurance,  high  approach  speeds  and  large  deck  space  requirements.  These 
characteristics  are  extremely  undesirable  traits  for  a  carrier-based  aircraft.  [Ref.  2] 


The  requirement  for  carrier  suitability,  therefore,  dictates  the  design  of  an  aircraft  with 
rounded,  moderately-swept  leading  edges.  When  such  a  planform  maneuvers  to  high 
angles  of  attack  the  vortices  produced  are  less  coherent  and  of  lower  energy.  Therefore, 
some  mechanism  is  needed  to  energize  or  induce  these  vortices  to  remain  coherent  at  high 
angles  of  attack,  without  the  added  penalty  of  poor  cruise  performance. 

B.  CANARD  CHARACTERISTICS 

Canards  are  separated  into  two  broad  categories:  long-coupled  or  close-coupled  (see 
Figure  2). 


Close-Coupled 


Figure  2.  Canard  Categories  [Ref.  2] 

Long-coupled  canards  are  of  the  type  used  primarily  as  a  control  surface  rather  than  as  a 
lifting  surface.  Examples  of  this  type  are  found  on  almost  all  missiles  and  on  some  aircraft, 
such  as  the  XB-70,  the  Concord  and  the  X-31  experimental  aircraft.  A  close-coupled 
canard  may  provide  a  significant  portion  of  the  aircraft  total  lift  in  addition  to  being  a 
control  surface.  The  aircraft  listed  at  the  beginning  of  the  introduction  are  examples  of 
close-coupled  designs. 

Extensive  research  has  been  performed  related  to  the  phenomenon  of  canard/wing 
interaction  [Ref.  2,3,6,71.  These  experimental  results  indicate  two  possible  mechanisms  by 


which  the  flowfield  from  the  canard  affects  the  flowfield  around  a  wing.  One  is  the 
impingement  of  the  canard  downwash  on  the  apex  of  the  main  wing;  the  other  is  the 
favorable  interference  between  the  canard  and  wing  leading-edge  vortices.  For  low  angles 
of  attack,  the  canard  produces  a  downwash  field  within  its  span  and  an  upwash  field 
outside  its  span.  The  result  is  a  nonuniform  angle  of  attack  on  the  main  wing.  The  inner 
(inboard  of  the  canard  span)  and  forward  portions  of  the  wing  have  a  lower  effective  angle 
of  attack  than  the  outer  and  rear  portions.  The  flow  over  the  wing  behind  the  canard  tends 
to  remain  attached  while  the  flow  outside  the  canard  or  at  the  rear  of  the  wing  tends  to 
separate.  The  probable  result  is  an  overall  loss  of  lift  on  the  main  wing  which  has  to  be 

compensated  for  by  the  increased  lift  of  the  canard.  In  fact,  studies  such  as  Reference  6 

o 
show  that  at  low  angles  of  attack  (<10  )  the  lift-curve  slope  of  a  model  with  or  without  a 

canard  is  identical.  The  second  mechanism  occurs  at  higher  angles  of  attack.  As  the  angle 
of  attack  is  increased,  a  strong  leading-edge  vortex  is  formed  on  the  canard.  As  this  vortex 
moves  over  the  wing  it  acts  to  energize  the  wing  leading-edge  vortex,  thereby  delaying 
Vortex-Breakdown  (VBD).  [Ref.  3,7] 
1.  Canard  Versus  Tail 

The  first  advantage  of  the  canard  over  the  horizontal  tail  design  arises  from  a 
difference  in  trim  requirements.  A  conventional  tail  balances  an  aircraft  in  flight  by 
producing  a  downward  lift  vector.  This  results  in  an  initial  decrease  in  lift  for  a  trim  to 
lower  speed.  A  canard,  on  the  other  hand,  produces  a  large  nose-up  pitching  moment 
which  must  be  balanced  by  a  positive  elevon  deflection  on  the  main  wing.  The  result  is 
increased  lift  for  the  canard  aircraft  due  simply  to  the  difference  in  the  trim  requirements. 
This  effect  is  most  dramatic  when  comparing  a  pure  delta-wing  aircraft  with  a  canard/delta- 
wing  combination.  Figure  3  shows  the  wing  loading,  approach  speed  and  lift  coefficient 
for  three  high-performance  aircraft.  The  Viggen  has  approximately  the  same  wing  loading 


as  the  F-106,  but  the  Viggen  approach  speed  is  34  knots  slower  than  that  of  the  F- 106  with 
a  70%  higher  lift  coefficient. [Ref.  2] 
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Figure  3.  Characteristics  of  the  F-4,F-106  and  Viggen  [Ref.  2] 

Additionally,  the  increased  lift  of  a  canard-configured  aircraft  cannot  be  completely 

accounted  for  simply  by  the  increased  lifting  area  of  the  canard.  Work  done  by  Er-El  [Ref. 

o 

3]  showed  an  increase  in  normal-force  coefficient  of  18%  at  22  angle  of  attack,  over  that 
of  a  corresponding  wing-alone  configuration.  The  canard  itself  added  only  9%  to  the  total 
lifting  area.  Likewise,  extensive  work  done  by  Lacey  [Ref.  2]  showed  similar  results. 
Figure  4(a)  shows  that  the  stall  of  a  wing  alone  and  of  a  wing/conventional-tail  combination 

o 
occurs  at  about  21  ,  whereas  there  is  no  indication  of  stall  for  the  canard-configured  aircraft 

o 
below  32  .  Additionally,  for  low  values  of  Ql  drag  of  the  all  three  configurations  was  the 

same,  see  Figure  4(b).  Therefore,  during  cruise  flight  the  canard/wing  configuration  would 


configuration.  Again,  the  increase  in  lift  of  the  canard/wing  configuration  at  high  angles  of 
attack  could  not  be  completely  accounted  for  by  the  additional  lifting  area  of  the  canard. 
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Figure  4a.    Lift  Coefficient  due  to  a  Canard  and  Horizontal  Tail  [Ref.  2] 
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Figure  4b.  Drag  Polar  for  a  Canard  and  Horizontal  Tail  [Ref.  2] 

C.  VORTEX  INTERACTION 

Er-El  [Ref.  3]  conducted  an  investigation  of  vortex  trajectory  and  breakdown  using 
canard/wing  configurations  at  angles  of  attacks  between  13°  and  30°.  Er-El  states  that  near 

the  wing  apex,  the  downwash  of  the  canard  forces  the  wing  vortex  down  closer  to  the 
wing,  while  further  downstream  the  vortex  trajectory  is  more  influenced  by  the  canard 
leading-edge  vortex  than  by  its  downwash.  That  is,  near  the  trailing  edge,  the  wing  vortex 
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is  displaced  upward  and  outward,  away  from  the  adverse  pressure  gradient.  Er-El  also 
states  that  a  close-coupled  canard  causes  a  delay  in  the  onset  of  the  breakdown  of  the  wing 
leading-edge  vortex,  which  originates  at  the  wing  apex.  This  delay  is  possibly  a  result  of 
the  wing  vortex  movement  away  from  the  adverse  pressure  gradient,  thereby  making  it 
more  stable. 

D.  PREVIOUS  TESTING 

At  the  Naval  Postgraduate  School,  Kersh  performed  force  measurements  on  a  model 
designed  in  accordance  with  Lacey's  guidelines.  [Ref.  6]  Kersh  investigated  the  effects  of 

O  0 

canard-enhanced  lift  at  five  angles  of  attack  between  10   and  48  .  Figure  5  [Ref.  6]  is  the 

lift  curve  that  resulted  from  that  study.  Note  that  as  mentioned  earlier  the  point  of  maximum 

lift  enhancement  occurred  at  22°,  the  point  of  first  stall  for  the  wing/body  configuration.  At 

o 
22    there  was  a  34%  increase  in  lift  for  the  canard/wing  configuration  over  that  for  the 

wing/body  configuration.  In  addition,  surface  flow  visualization  indicated  the  dramatic 
effect  the  canard  vortex  had  on  the  main  wing  flowfield.  However,  while  the  flow 
visualization  gave  outstanding  insight  into  the  nature  of  the  flow  on  the  surface  of  the  wing, 
it  did  not  say  much  about  what  occurred  off  the  wing  surface.  That  flowfield  study  was  the 
major  objective  of  this  study. 
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Figure  5.  Lift  Coefficient  of  Wing/body  and  Maximum  Lift  of  Canard/Wing 

[Ref.  6] 
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III.  EXPERIMENT  AND  PROCEDURES 

A.   PURPOSE 

This  study  was  conducted  to  investigate  the  effects  of  vortex  interference  on  the 
flowfield  around  the  main  wing  of  a  close-coupled  canard  and  wing  combination.  A  nulling 
five-hole  pressure  probe  was  used  to  make  flowfield  pressure  measurements  in  the  three 
crossplanar  grid  locations  shown  in  Figure  6.  These  pressure  measurements  were  then 
used  to  determine  local  flow  velocities,  pressures,  and  pitch  and  yaw  angles  for  use  in 
mapping  the  position  and  interaction  of  canard  and  wing  vortices. 
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Figure  6.  Model  Overview 
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B.    APPARATUS 

The  primary  equipment  used  in  this  study  were  a  wind  tunnel,  a  close-coupled  canard 
model,  a  three-dimensional  traversing  mechanism,  a  rotary  pressure  transducer,  a  data 
acquisition  system,  and  a  nulling  five-hole  pressure  probe. 

1.  Wind  Tunnel 

The  tunnel  was  of  the  close-circuit,  single-return  type.  It  measured  64  feet  in  length 
and  between  21.5  and  25.5  feet  in  width  with  a  test  section  cross-sectional  area  of  8.75 
square  feet.  The  tunnel  was  powered  by  a  lOOhp  electric  motor  coupled  to  a  three-bladed 
variable  pitch  fan  via  a  four-speed  transmission.  Figure  7  is  a  schematic  of  the  tunnel. 
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Figure  7.  Naval  Postgraduate  School  Wind  Tunnel  [Ref.  8] 
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Tunnel  turbulence  was  minimized  by  the  use  of  stator  blades,  turning  vanes  and 
turbulence  screens.  Eight  stator  blades  were  used  to  remove  swirl  imparted  to  the  flow  by 
the  fan  blades.  Turning  vanes  were  located  at  each  corner  to  turn  the  flow  while 
maintaining  its  kinetic  energy  and  turbulence  screens  were  placed  in  the  settling  chamber. 

The  test  section  measured  45  inches  by  28.5  inches  with  a  cross-sectional  area 
about  one  tenth  that  of  the  settling  chamber.  The  test  section  was  rectangular  with  slightly 
divergent  walls  to  account  for  boundary  layer  growth.  Breather  slots  at  the  far  end  of  the 
test  section  were  used  to  keep  the  test  section  at  approximately  atmospheric  pressure.  Due 
to  the  configuration  of  the  pressure  measuring  equipment,  ambient  atmospheric  pressure 
was  used  as  the  reference  static  pressure  rather  than  test  section  static  pressure.  [Ref.  8] 

Tunnel  velocity  was  set  by  reference  to  a  water  manometer  that  measured  the 
difference  in  static  pressure  between  the  settling  chamber  and  the  test  section.  The  static 
pressures  in  the  settling  chamber  and  at  the  test  section  entrance  were  determined  by 
reference  to  four  pressure  ports,  one  on  each  wall.  The  pressures  from  these  ports  were 
averaged  by  a  common  manifold  prior  to  the  manometer.  The  manometer  measured  the 
pressure  difference  in  centimeters  of  water  and  the  test  section  velocity  was  then  determined 
by  equation  (1). 

V  =  {(2*2.0475*PcmH2O)/p*K)}  1/2  (1) 

where: 

V  Test  Section  Velocity  (ft/sec) 

2.0475  Converts  From  Centimeters  of  Water  to  Lbf  /  ft2 

Pcm^O  Manometer  Reading  (cm  of  water) 

p  Density  of  Air  (slugs/ft^) 

K  Tunnel  Calibration  Factor  (PcmH20/q) 
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The  tunnel  calibration  factor,  the  ratio  of  the  static-pressure  ratio  to  the  test  section 
dynamic  pressure,  was  determined  to  be  0.925. 
2.   Model  Design 

The  model  used  in  this  study  was  originally  used  by  Kersh  [Ref.  6]  and  was 
designed  in  accordance  with  the  earlier  work  of  Lacey  [Ref.  2].  During  Lacey's  studies,  a 
canard-area-to-wing-area  ratio  of  0.20  as  referenced  to  the  model  centerline  of  the  fuselage 
was  used.  Because  it  was  desired  to  accommodate  an  internal  mechanism  for  a  movable 
canard,  the  fuselage  of  this  model  was  larger  than  that  of  Lacey's  model.  Therefore,  this 
model  was  designed  with  an  exposed  canard-area-to-wing-area  (wing  referenced  to  the 
centerline  of  the  fuselage)  ratio  of  0.20.  If  the  area  ratio  had  not  been  determined  this  way, 
a  canard  exposed  area  much  smaller  than  Lacey's  would  have  resulted.  If  the  area  ratio  of 
Lacey's  model  had  been  determined  by  this  same  method,  a  ratio  of  0.13  would  have 
resulted.  Consequently,  an  exposed  canard-area-to-wing-area  ratio  of  0.20  compared 
favorably  with  that  of  Lacey's  model. 

Lacey's  work  also  showed  that  the  relative  position  of  the  canard  to  the  wing  was 
crucial  for  constructive  canard/wing  interference.  According  to  Reference  2,  the 
longitudinal  position  of  the  canard  (X),  non-dimensionalized  by  the  mean  aerodynamic 
chord  (MAC)  of  the  wing  referenced  to  the  fuselage  centerline  (X/CmacX  should  not  be 
greater  than  X/Cmac=l-5.  The  longitudinal  distance  was  measured  from  the  40%  chord 
location  of  the  exposed  canard  root  to  the  quarter  chord  of  the  wing  MAC.  Also,  there 
should  be  no  overlap  of  the  main  wing  and  canard.  As  a  result,  the  canard  was  positioned 
at  X/Cmac=l-2.  The  canard  vertical  position  was  Z/Cmac=0-2;  which  was  chosen  such 
that  it  would  match  Lacey's  model. 

The  wing  and  canard  were  both  an  NACA  64A008  airfoil  section,  which  was  the 

o 

same  section  used  by  Lacey.  The  canard  had  a  leading-edge  sweep  of  60  ,  a  straight-taper 

o 

of  0. 1  and  an  aspect  ratio  of  2.0.  The  wing  had  a  leading-edge  sweep  of  50  ,  a  straight- 
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taper  of  0.15  and  an  aspect  ratio  of  3.0.  These  planforms  were  chosen  to  ensure  that  strong 
leading-edge  vortices  would  be  generated.  Unlike  the  sharp  leading-edge  model  used  in  Er- 
El's  studies  [Ref.  3],  both  the  canard  and  wing  had  rounded  leading  edges.  No  attempt  was 
made  to  trip  the  boundary  layer.  Planform  geometry  was  derived  from  equations  (2),  (3) 
and  (4).  [Ref.  6] 

AR  =  2{b/Cr(l+>.)}  (2) 

AR  =  b2/S  (3) 

MAC  =  2/3  { Cr  +  Q  -  [CrCt/(Cr  +  Q)] }  (4) 


Where: 

AR 

Aspect  Ratio 

b 

Span 

Q 

Length  of  Root  Chord 

Q 

Length  of  Tip  Chord 

X 

Taper  Ratio  (Ct/Cr) 

s 

Area  of  Wing  or  Canard 

MAC     Wing  Mean  Aerodynamic  Chord 

Figure  8  is  a  schematic  of  the  model  and  contains  a  listing  of  its  geometric  data. 
Reference  6  contains  a  thorough  description  of  the  model's  design.  Because  the  five-hole 
pressure  probe  was  not  long  enough  to  reach  the  model  mounted  on  the  tunnel  floor,  a 
platform  was  constructed  to  elevate  the  model  (see  Figure  9).  The  effect  on  the  flowfield 
should  be  minimal  for  the  half-model,  as  the  raised  platform  acted  in  a  similar  manner  to  the 
reflection  plane. 
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Figure  8.  Model  Geometric  Data 
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Figure  9(a).  Side  View  of  the  Model  Mounted  in  the  Tunnel 
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Figure  9(b).  Front  View  of  the  Model  Mounted  in  the  Tunnel 
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3.  Velmex  8300  Three-Dimensional  Traverser  System 

The  Velmex  system,  shown  in  Figure  10,  was  used  to  accurately  position  the  five- 
hole  pressure  probe  in  the  test  section.  The  system  consisted  of  the  traverser  mechanism, 
Figure  10(b),  and  the  motor  controller,  Figure  10(a).  The  traverser  mechanism  was 
composed  of  leadscrews,  slides  and  motor/jackscrew  assemblies.  Each  motor  was  a  200- 
step-per-revolution,  10-amp  stepping  motor  with  a  maximum  velocity  of  3000  steps  per 
second.  The  motor  step  size  was  1/200  of  a  revolution  which  equated  to  0.000125  inch. 
The  three  motor  assemblies  received  their  commands  from  the  controller  which  was  capable 
of  interpreting  signals  from  either  a  parent  computer  or  manual  inputs  at  the  controller. 
During  this  study  the  controller  received  its  commands  from  a  parent  PC/AT. [Ref.  9] 

4.  Rotary  Pressure  Transducer  and  Data  Acquisition  System 

A  48-port  rotary  pressure  transducer  (commonly  referred  to  as  a  Scanivalve®)  read 
each  port  of  the  pressure  probe  twice  per  measurement  point.  The  data  acquisition  system, 
shown  in  Figure  11,  consisted  of  the  required  hardware  and  software  to  allow  a  PC/AT  to 
acquire  the  pressure  data.  The  Scanivalve®  put  out  a  7-bit  binary  coded  decimal  signal  that 
corresponded  to  the  port  currently  being  monitored.  This  arrangement  allowed  remote, 
electronic  monitoring  and  control  of  the  Scanivalve®  assembly.  [Ref.  10]  After  the 
pressure  signal,  in  the  form  of  an  analog  voltage,  left  the  Scanivalve®  it  was  passed  to  the 
Relay  Multiplexer,  then  to  the  Digital  Multimeter  (DMM).  The  DMM  converted  the  analog 
voltage  signals  to  a  digital  form  which  were  then  sent  to  the  computer  for  further 
processing. 
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Figure  10(a).     Velmex  8300  Traverser  -  Motor  Controller 


Figure  10(b).     Velmex  8300  Traverser  Assembly 
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Figure  11.  Data  Acquisition  Hardware  [Ref  10] 
5.  Five-Hole  Pressure  Probe 

The  three-dimensional  five-hole  pressure  probe,  shown  in  Figure  12,  was  made  of 
corrosion-resistant,  non-magnetic,  stainless  steel.  It  was  0.125  inches  in  diameter  and  24 
inches  in  total  length.  At  the  tip  of  the  probe  were  five  measuring  holes.  The  center  hole  of 
the  probe  (P^)  measured  total  pressure,  the  two  lateral  holes  (P2»P3)  were  used  to  measure 
yaw  angle  and  static  pressure,  and  the  top  and  bottom  holes  (P^Ps)  were  used  for  pitch 
information.  [Ref.  11] 


Block  Size 

1.75"  x  0.5" 


Figure  12.  The  Five-Hole  Pressure  Probe  [Ref.   11] 
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C.  EXPERIMENT  SOFTWARE 

The  software  used  in  this  study  was  developed  by  Lung  during  his  work  with  the  five- 
hole  pressure  probe  [Ref.  10].  Some  modifications  were  made  to  accommodate  this  work, 
but  the  flow  of  data  remained  essentially  the  same.  Figure  13  is  a  schematic  flowchart  of 
the  various  programs  and  data  files.  The  following  sections  discuss  the  individual 
programs. 

1.   PPROBE 

PPROBE  (Appendix  A),  written  in  BASICA  by  Kindelspire  [Ref.  12],  and 
modified  by  Lung  [Ref.  10],  was  used  to  control  the  data  acquisition  process.  Figure  14  is 
the  program  flow  chart  for  PPROBE.  The  two  primary  responsibilities  of  PPROBE  were 
the  control  of  the  traverser  and  the  recording  of  the  data. 

First,  PPROBE  was  used  to  manually  position  the  probe  to  the  proper  starting  point 
in  the  tunnel,  then  the  dimensions  of  the  grid  to  be  measured  and  the  desired  step  size  were 
entered.  PPROBE  calculated  the  total  number  of  points  to  be  measured,  created  the  required 
number  of  data  files  and  initialized  the  data  acquisition  equipment.  The  operator  then 
physically  nulled  the  probe  and  entered  the  observed  yaw  angle  into  PPROBE.  PPROBE 
rotated  the  pressure  transducer  to  the  appropriate  port  and  began  measuring  the  analog 
voltage  which  corresponded  to  the  pressure  at  that  port.  Finally,  the  data  were  displayed  to 
the  operator  and  he  was  given  a  chance  to  remeasure  the  data  point.  If  the  data  were 
acceptable  PPROBE  moved  the  traverser  to  the  next  measuring  point.  At  the  completion  of 
each  column  in  the  grid,  the  data  were  placed  into  a  file  to  be  used  in  the  program  PVA. 
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Figure  13.  Data  Flow  Path 
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Figure  14.  PPROBE  Flow  Chart 
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2.  CALP 

CALP  (Appendix  B),  was  written  in  BASICA  and  used  to  calibrate  the  pressure 
transducer  prior  to  each  data  collection  period.  The  calibration  manometer  is  shown  in 
Figure  15.  Six  pressures  were  recorded  using  the  calibration  manometer  and  the  pressure 
transducer  in  the  Scanivalve®  to  provide  a  known  calibration  curve  for  the  experiment.  A 
straight  line  was  fit  to  the  data,  whose  slope  and  intercept  were  input  to  the  PVA  program. 


Figure  15.  Calibration  Manometer 
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3.  PVA  (Pressure  Velocity  Angle) 

PVA  (Appendix  C)  was  a  modified  version  of  CONVERT,  a  Fortran  program 
which  combined  the  output  of  PPROBE  and  CALP  to  calculate  the  flowfield  properties 
[Ref  10].  PVA  used  the  pressure  data,  in  the  form  of  analog  voltages,  and  yaw  data  from 
PPROBE,  and  combined  them  with  the  curve  fit  derived  from  the  CALP  output  to 
determine  the  desired  flowfield  characteristics.  See  Figure  16  for  the  PVA  flow  chart. 
Specifically,  PVA  determined: 

-freestream  velocity  (V) 

-crossflow  components  of  freestream  velocity  (Vx  and  Vy) 

-total  pressure  and  total  pressure  coefficients 
-static  pressure  and  static  pressure  coefficients 
-flow  yaw  and  pitch  angles 
The  pressure  coefficients  were  determined  as  follows: 

Cpr={PTl-Pamb-Ql}/Ql  (5) 

Cps  =  {Psl-Pamb}/Ql  (6) 

Where: 

Pri    =  Local  Dynamic  Pressure 
Psl     =  Local  Static  Pressure 
Pamb  =  Room  Ambient  Pressure 
Ql      =  Freestream  (Tunnel)  Dynamic  Pressure 
Due  to  the  configuration  of  the  pressure-measuring  equipment,  room  ambient  pressure  was 
used  as  the  reference  static  pressure  in  the  pressure  coefficient  calculations.  The  output  of 
PVA  was  put  in  an  ASCII  file  for  use  by  commercially  available  graphic  software  for 
plotting  the  desired  output. 
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PVA  Flow  Chart  (1) 
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Figure  16(a).  PVA  Program  Flow  Chart 
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From  PVA  (1) 
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Figure  16(b).  PVA  Program  Flow  Chart 
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D.  EXPERIMENTAL  CONDITIONS  AND  PROCEDURES 

To  correlate  these  results  with  the  results  of  Kersh,  the  experimental  conditions  were 
kept  as  constant  as  possible.  A  reference  measured  pressure  difference  of  17  cm  of  H2O 
was  set,  resulting  in  a  test-section  dynamic  pressure  of  37.63  lbf/fn  and  an  average 
velocity  of  177  ft/s.  The  Reynolds  number  was  Re=9.5xlCr  based  on  the  main  wing  mean 

aerodynamic  chord. 

To  account  for  possible  tunnel  crossflow,  probe  installation  effects  and  flowfield 
disturbances  caused  by  the  model  platform,  an  initial  run  was  made  without  the  model 
installed,  to  established  baseline  yaw  and  pitch  values.  This  measurement  was  again  run  at 
the  end  of  testing  to  ensure  the  baseline  values  had  not  changed  appreciably. 

Measurements  were  made  at  three  locations,  numbered  as  shown  in  Figure  6.  Grids  1 
and  2  were  mapped  at  1/4"  intervals,  while  grid  3  was  mapped  at  1/2"  intervals.  Note  that 
the  measuring  plane  of  the  grids  were  perpendicular  to  the  freestream  direction,  rather  than 
to  the  wing  chord. 

Grid  1  was  4"  by  5"  and  was  set  just  behind  the  canard  trailing  edge  in  order  to  capture 
the  canard  leading-edge  vortex.  Grid  number  2  was  5.5"  by  6"  and  was  placed  at  the  mid- 
point of  the  wing  to  capture  the  effect  of  the  canard  leading-edge  vortex  on  the  flowfield  at 
that  point.  Finally,  since  flow  visualization,  completed  earlier,  had  shown  the  flowfield 
beginning  to  lose  its  coherency  near  the  trailing  edge  of  the  wing,  grid  number  3  was 
placed  just  behind  the  trailing  edge  of  the  wing.  Grid  3  measured  7.5"  by  6";  because  its 
size  would  have  necessitated  an  excessive  amount  of  tunnel  time  to  complete  at  1/4" 
intervals,  it  was  mapped  at  1/2"  intervals. 
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IV.   RESULTS 

The  following  sections  discuss  the  results  of  this  study.  In  each  section  there  will  be 
three  plots.  The  first  will  be  a  plot  of  flowfield  velocity  vectors,  then  a  plot  of  the  velocity 
streamlines,  followed  finally  by  a  plot  of  total-pressure-coefficient  contours.  Flow 
visualization  results  are  also  included  to  correlate  with  the  measured  data.  All  the  plots  are 
scaled  in  inches,  with  "X"  plotted  horizontally,  transversely  to  the  tunnel  and  "Y"  plotted 
vertically  to  the  tunnel.  Locations  on  the  plots  will  be  identified  by  listing  the  coordinates  in 
the  standard  (X,Y)  convention.  Also,  the  plotted  data  are  shown  with  cutaway  views  of  the 
model  geometry  to  establish  the  relative  position  of  the  measuring  plane  to  the  model.  At 
grid  1  the  actual  measuring  plane  was  just  beyond  the  canard  trailing  edge,  but  the  canard 
cross-section  at  the  trailing  edge  will  be  shown  for  reference.  Likewise,  at  station  three  the 
grid  was  located  after  the  main  wing  trailing  edge,  but  the  wing  cross  section  3  is  shown  to 
indicate  the  relative  position  of  the  vortices  and  main  wing  (see  Figure  17). 


Figure  17.  Model  Side  view 
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Because  of  the  dynamic  nature  of  flowfields  surrounding  aerodynamic  surfaces  at  high 
angles  of  attack,  there  were  regions  in  which  the  capabilities  of  the  probe  were  exceeded.  In 
particular,  the  probe  was  unable  to  make  accurate  measurements  downstream  of  stalled 
aerodynamic  surfaces  or  near  the  core  of  a  vortex.  Precisely,  in  these  regions  the  probe 
could  not  be  nulled  -  that  is,  P2  and  P3  could  not  be  made  equal.  The  main  effect  of  this 
probe  limitation  was  that  the  exact  positions  of  the  vortex  cores  could  not  be  determined. 
During  data  reduction,  to  account  for  the  regions  of  the  flowfield  where  the  probe  could  not 
be  nulled,  a  value  of  zero  velocity  was  assigned  and  a  limiting  value  of  total  pressure  drop 
was  established  based  on  observed  trends  in  the  area  of  the  flowfield  where  accurate  data 
were  available.  Therefore,  the  velocity  vector  plots  and  velocity  streamline  plots  have  areas 
where  there  are  no  data  present;  these  areas  are  the  sections  of  the  flowfield  where  the 
probe  capabilities  were  exceeded.  Likewise,  on  the  pressure  contour  plots  an  area  of 
uniform  total  pressure  loss  indicates  the  section  in  which  the  data  are  inaccurate.  The  isobar 
along  which  the  pressure  data  become  unreliable  will  be  delineated  in  the  section  which 
refers  to  that  plot. 

An  additional  problem  of  the  probe's  lack  of  capabilities  arises  from  the  fact  that  both 
vortical  flow  and  the  reversed  flow  from  a  stalled  aerodynamic  surface  appear  as  blank 
areas  in  the  data  or  as  areas  of  total  pressure  loss.  Since  there  is  insufficient  data  to  indicate 
a  difference  between  the  two  flow  states,  flow  visualization  results  have  been  included  to 
discern  between  regions  of  separated/reversed  flow  and  those  of  separated  coherent  vortex 
flow. 

A.   WING-ALONE 

With  the  canard  off,  the  flow  over  the  main  wing  was  characterized  by  large  regions  of 
reversed  flow,  indicative  of  a  stalled  wing.  In  general,  as  can  be  seen  from  Figure  18,  there 
was  a  large  amount  of  oil  over  the  majority  of  the  wing  indicating  wing  stall.  Particularly, 
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the  oil  on  the  inboard  portion  of  the  wing  was  mostly  stagnate  or  flowing  spanwise;  this 
response  indicated  that  this  section  of  the  wing  had  stalled  and  thereby  created  an  area  of 
large  total  pressure  loss.  Along  the  leading  edge  of  the  wing  there  appeared  to  be  a  weak, 
poorly  organized,  separated  vortex,  which  apparently  lost  its  coherency  over  the  outboard 
section.  This  separated  vortex  was  indicated  in  Figure  18  by  the  heavy  line  of  oil  which 
started  at  the  wing  apex  and  then  seemingly  burst  at  the  wing  mid-point.  Therefore,  it 
would  be  expected  that  at  each  measuring  grid  there  would  be  evident  a  general  tendency 
for  the  flow  to  move  spanwise,  in  a  root-to-tip  fashion,  with  a  large  area  of  total  pressure 
loss.  In  this  case,  the  loss  in  total  pressure  would  be  due  to  chaotic,  disorganized, 
separated  reversed  flow  rather  than  coherent,  vortical  flow. 


Figure  18.  Wing  Alone  Set  at  22 
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1.  Mid-Point  of  the  Main  Wing  (Grid  2) 

Figure  19(a)  is  a  plot  of  the  "X"  and  "Y"  components  of  the  local  velocities  at  grid 
2.  The  blank  area  inside  the  heavy  line  indicates  the  area  where  the  data  were  unreliable  due 
to  flow  separation.  The  crossflow  component  of  freestream  velocity  increased  as  the  probe 
was  traversed  from  the  wing  tip  to  the  root.  The  crossplane  velocity  streamline  plot  of 
Figure  19(b)  shows  a  vortical  pattern  which  compares  favorably  with  the  flow  visualization 
of  Figure  18.  Crossplane  velocities  are  vertical  (in  the  model  reference  sense)  near  the 
leading  edge  or  tip,  but  horizontal  at  26%  mean  aerodynamic  chord  (MAC)  above  the  wing. 
Figure  20  is  a  contour  plot  of  total-pressure-coefficient  loss.  This  plot  also  shows  an  area 
of  total  pressure  loss  increasing  in  size  as  the  probe  was  moved  inboard  toward  the  wing 
root.  At  the  outboard  section  of  the  wing  it  is  possible  that  the  cause  of  this  loss  was  the 
disintegrating  vortex  that  appeared  to  form  there;  however,  at  the  inboard  section,  the  loss 
was  certainly  caused  by  reversed  flow.  Pressures  inside  the  -1.24  isobar  are  unreliable. 
The  separated  region  extends  up  to  29%  MAC. 

2.  Trailing  Edge  of  the  Main  Wing  (Grid  3) 

At  grid  3,  which  was  located  just  past  the  trailing  edge  of  the  wing,  Figures  21(a) 
and  (b)  clearly  indicated  a  strong  spanwise  flow  pattern  from  the  wing  root  to  the  tip  near 
the  fuselage.  This  flow  pattern  is  believed  to  be  caused  by  entrainment  due  to  the  separated 
reversed  flow  commonly  associated  with  wing  stall,  rather  than  by  a  coherent  vortex.  The 
separation  now  extends  about  57%  of  the  MAC  above  the  wing  plane  at  the  trailing  edge. 
The  total-pressure-coefficient  gradients  in  Figure  22  show  that  the  region  of  total  pressure 
loss,  due  to  the  stalled  reversed  flow,  had  grow  substantially  as  the  flow  moved 
downstream  from  grid  2.  Pressures  inside  the  -1.20  isobar  are  considered  to  be  unreliable. 
Though  the  region  of  uncertain  data  is  large,  the  results  seem  to  indicate  the  extent  to  which 
the  flow  separation  affects  the  wing  flowfield,  just  as  the  oil  flow  visualization  seemed  to 
indicate  the  same  on  the  wing  surface. 


32 


Y  (inches) 


o 

c 


X 


■a 

u 
as 
c 
« 
U 

o 
Z 

OX) 

c 


4» 


o 
o. 


u 
O 

■4—1 
O 


s5 


L. 

3 

ex) 


33 


Y  (inches) 


e/3 

(J 

c 


X 


T3 
u 

« 

C 

r? 

u 

o 

z 

OX) 

c 


o 


c 


t- 
(73 


3 
OX) 


34 


<*)a>mo<oc\ioo''»'a>m>-i--(MoOT 

T-(tj«iosoiq(\iniflNaoT-n 

000000»-T-i-r-»-*-c\JC\JC\J 


■a 

L. 

a 
c 
« 

u 

o 

0X1 

c 


4> 


o 

Q. 

1/5 

u 
3 
O 

C 

o 
U 


c 
"a 

<x> 

O 

U 

U 

3 

Vi 

ft 

4> 
U 


o 
H 

© 

4» 
U 
3 
OX) 


35 


Y  (inches) 


o 


q 
in 


o 


q 
06 


q 
cvi 


q 
d 


c 


X 


q 
cd 


o 


q 

CO 


q 
c\i 


q 
d 


1      1 

| 

I 

1 

1 

I 

1 

1 

I 

\ 

\ 

V 

X 

V 

X 

^ 

* 

V 

^ 

^ 

X 

X 

X 

X 

X 

X 

\ 

^ 

\ 

\ 

\ 

* 

*            X 

X 

X 

x 

X, 

^ 

V 

\ 

\ 

\ 

\ 

Vl 

»                    X 

•% 

^ 

^ 

•x 

^ 

\ 

\ 

X 

X 

X 

X 

- 

»               •* 

X 

X 

x. 

X 

-*. 

X 

X 

\ 

\ 

* 

* 

- 

^ 

^ 

^ 

^^. 

^ 

^ 

"^ 

*        ~^s 

*         , 

— 

i        , 

*        \ 

'         I 

- 

'        \ 

- 

*        1 

" 

/ 

I 

■o 
u 
<a 
c 
a 

U 

o 
Z 

ox 


o> 


■o 
■ 

C 

"3 

u 
H 


en 

u 

o 


o 
"a* 


«J 


j- 
3 
0£ 


36 


Y  (inches) 


W5 

o 

c 


X 


■o 
u 
« 
c 
« 

U 

o 
z 

c 


4» 
OX) 

■a 
■ 

OX) 

c 


en 
H 


I 

« 

U 

■*■> 
CO 


u 

3 
OX) 


37 


T3 

s_ 

C 

U 

o 
Z 

C 


0> 


4> 

M) 

w 

OH 
C 

"« 

u 

H 

u 
3 
O 
■*-* 

e 
o 
U 


c 

o» 
o 
U 

o> 

u 
3 

1/3 
<M 

o> 

u 


O 

H 


u 
3 


38 


B.  WING/CANARD  COMBINATION 

The  addition  of  the  canard  to  the  model  caused  large-scale  reattachment  of  the  flow  over 
the  wing.  Figure  23  shows  a  well-defined  secondary  separation  line  along  the  leading 
edges  of  both  the  canard  and  wing.  These  separation  lines  are  characteristic  of  strong, 
highly  coherent  leading-edge  vortices.  Note  in  Figure  23  that  despite  the  high  angle  of 
attack  of  the  canard  (29°),  the  flow  appeared  to  remain  completely  attached  (not  separated 
and  reversed)  and  that  the  large  leading-edge  sweep  (60°)  produced  a  strong  vortex  even 
though  the  leading  edge  was  rounded.  In  reference  to  Figure  23,  the  wing  can  be  seen  to  be 
fairly  clean  of  oil  with  the  exception  of  a  wedge  which  started  just  aft  of  the  apex  and  grew 
as  it  proceeded  downstream  parallel  to  the  fuselage.  Also,  there  is  a  line  of  oil  along  the 
trailing  edge  of  the  wing  which  indicated  that  the  flow  had  begun  to  separate  there, 
particularly  at  the  point  were  the  trailing-edge  line  joined  the  line  which  had  started  near  the 
apex  and  moved  parallel  to  the  fuselage.  Additionally,  looking  at  the  top  surface  of  the 
fuselage  where  the  canard  and  wing  join  the  body,  one  can  see  there  is  a  section  of  the 
fuselage  free  of  oil.  It  is  thought  that  this  result  was  due  to  a  secondary  vortex  forming  near 
the  surface  of  the  fuselage.  The  exact  origin  of  this  secondary  fuselage  vortex,  however,  is 
not  clear.  It  could  have  been  the  result  of  either  the  forebody  or  of  the  leading-edge  vortices 
formed  on  the  canard  and  wing,  or  possibly  some  combination  of  the  two  effects. 

The  discussions  that  follow  will  be  in  the  same  format  used  in  the  preceeding  wing- 
alone  discussion.  The  path  of  the  canard  and  wing  vortices  will  be  analyzed  and  the  effects 
of  the  one  on  the  other  noted. 
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0  0 

Figure  23.  Model  set  at  22  ,  Canard  set  at  7     Incidence 


1.  Trailing  Edge  of  the  Canard  (Grid  1) 

The  flow  at  this  grid,  which  was  located  just  aft  of  the  canard  trailing  edge,  was 
dominated  by  a  strong  vortex  generated  at  the  canard  leading  edge.  Figure  24(a)  shows  that 
for  a  large  portion  of  the  measuring  grid,  the  probe  was  unable  to  accurately  record  flow 
field  pressures.  However,  Figures  24(a)  and  (b)  do  indicate  the  presence  of  very  strong 
vortical  flow  in  a  clockwise  direction.  Also,  Figures  24(a)  and  (b)  indicate  the  presence  of  a 
secondary  vortex  forming  toward  the  fuselage  surface  near  the  canard/body  juncture.  The 
secondary  vortex  was  rotating  in  a  counter-clockwise  direction,  which  indicated  that  it  had 
probably  formed  as  a  result  of  the  canard  leading-edge  vortex  in  much  the  same  manner  as 
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a  secondary  vortex  forms  near  a  wing  leading-edge  vortex.  The  mechanism  for  such  a 
vortex  would  be  the  same  as  that  depicted  in  Figure  1,  but  would  be  a  result  of  an 
interaction  of  the  wing  and  fuselage  rather  than  the  wing  and  itself.  Such  a  phenomenon 
had  not  been  noted  in  previous  testing  because  the  models  used  did  not  have  a  substantial 
fuselage.  If  this  secondary  vortex  had  resulted  from  the  model  forebody,  it  would  have 
been  expected  to  be  rotating  clockwise. 

Figure  25  shows  the  steep  pressure  gradients  which  resulted  from  this  vortex.  The 
data  inside  the  -1.26  isobar  are  unreliable,  but  clearly,  the  trend  was  for  a  tight,  coherent 
core  vortex  standing  approximately  1.5  to  2.5  inches  off  the  canard  surface  at  the  trailing 
edge.  The  presence  of  the  secondary  vortex  near  the  canard  root  is  better  indicated  in  the 
plot  of  total-pressure-coefficient-loss  contours.  Notice  that  a  freestream  value  of  total 
pressure  extends  all  the  way  to  the  canard  surface,  which  agrees  with  the  flow  visualization 
of  Figure  23. 
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2.  Mid-Point  of  the  Main  Wing  (Grid  2) 

A  comparison  of  Figures  26  and  19  shows  a  dramatic  difference  in  the  flowfield 
due  to  the  addition  of  the  canard  to  the  model.  Figure  26(a)  shows  two  regions  of 
unobtainable  data,  one  at  approximately  (5,2)  and  the  other  very  near  the  wing  leading  edge 
at  about  (0.5,3.5).  The  area  located  near  (5,2)  was  the  canard  leading-edge  vortex  which 
had  moved  downstream  from  grid  1.  The  area  near  (0.5,3.5)  could  have  been  the  result  of 
two  possible  phenomena.  One  is  that  the  large  region  of  chaotic  flow  seen  in  Figure  19  had 
been  forced  into  a  small  "pocket"  of  separated  reversed  flow.  The  second  is  that  the  wing 
leading-edge  vortex  had  reformed  due  to  the  presence  of  the  canard  vortex.  The  flow 
visualization  in  Figure  27  indicates  the  presence  of  a  wing  leading-edge  vortex. 
Furthermore,  the  pattern  formed  by  the  streamlines  in  Figure  26(b)  indicates  that  rather  than 
the  flow  having  merely  been  forced  into  a  "separation-pocket",  the  leading-edge  vortex  was 
actually  energized  by  the  canard  vortex  and  had  reformed.  The  flow  visualization  suggests 
that  a  strong  vortex  has  formed  along  the  wing  surface. 

The  energizing  of  the  wing  leading-edge  vortex  in  this  manner  and  the  reattached 
flow  over  the  inboard  wing  section  would  lead  to  a  level  of  enhanced  lift  which  would  be 
greater  than  that  obtained  by  simply  adding  the  lifting  surface  of  the  canard.  In  fact,  as  was 
mentioned  earlier,  this  was  found  to  be  exactly  the  case  in  previous  studies  [Ref.  2,3].  The 
mechanism  is  clearly  demonstrated  here.  Also  in  Figure  26  can  be  seen  a  significant 
crossflow  component  between  the  canard  vortex  and  the  wing  leading-edge  vortex.  In  the 
vicinity  of  (2,1)  the  crossflow  is  as  much  as  65%  of  the  total  freestream  velocity  and 
approaches  the  wing  perpendicular  to  its  surface.  The  formation  of  oil  parallel  to  the 
fuselage  in  Figure  27  could  have  been  due,  at  least  in  part,  to  the  stagnation  of  this  large 
crossflow  component  on  the  wing  surface. 

The  total-pressure-coefficient  contours  of  Figure  28  show  the  presence  of  a  third 
vortex  not  completely  captured  by  the  velocity-vector  or  velocity-streamline  plots.  Centered 
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at  (2.1,0.6),  this  vortex  is  believed  to  be  the  secondary  vortex  that  had  formed  on  the 
fuselage  surface  at  grid  1.  Most  importantly,  the  contours  of  Figure  28  show  that 
freestream  total  pressure  extended  to  the  wing  surface  over  the  inboard  section.  Clear 
evidence  is  provided  that  the  flow  had  reattached  over  the  inboard  portion  of  the  wing, 
enhancing  significantly  the  lifting  ability  of  the  previously  stalled  wing. 
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Figure  27.  Wing  at  22°  with  Canard  Deflected  7° 
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3.  Trailing  Edge  of  the  Main  Wing  (Grid  3) 

While  the  velocity  vectors  and  streamlines  of  Figure  29  indicate  a  similar  pattern  to 
that  of  grid  2,  the  flow,  in  general,  was  less  coherent  and  slightly  more  diffuse  than  at  grid 
2. 

Figure  29(a)  shows  that  there  were  two  areas  devoid  of  obtainable  data.  One  of 
these  areas  was  the  result  of  the  canard  vortex  which  had  moved  downstream  to  a  point 
near  the  (5.5,1.5)  position  and  the  other  was  caused  by  the  wing  leading-edge  vortex 
which  had  reformed  as  a  result  of  the  presence  of  the  canard.  The  streamlines  in  Figure 
29(b)  show  that  the  flow  at  these  locations  was  strongly  vortical.  In  addition,  the  flow 
visualization  of  Figure  27  shows  the  presence  of  a  secondary  separation  line  consistent 
with  the  formation  of  a  leading-edge  vortex.  Figure  29  also  indicated  that  there  still  existed 
a  strong  crossflow  component  moving  from  the  canard  vortex  to  the  wing  vortex.  The 
magnitude  of  this  crossflow  reached  66%  of  the  freestream  as  at  grid  2;  however,  the 
freestream  velocity  on  the  wing  surface  of  grid  3,  near  (0,2.5),  had  slowed  to  74%  of  that 
at  the  same  position  on  grid  2.  This  was  substantiated  by  the  growth  in  the  width  of  the  oil 
line  parallel  to  the  fuselage.  A  comparison  of  Figures  18  and  21  with  29  shows  that  the 
energy  introduced  to  the  flow  by  the  canard  vortex  moved  the  burst  point  of  the  wing 
leading-edge  vortex  past  the  trailing  edge  of  the  wing. 

The  pressure  contours  in  Figure  30  show  the  presence  of  the  two  vortices  already 
mentioned  and  a  third  area  of  pressure  loss.  This  third  area  was  most  likely  a  combination 
of  the  remnants  of  the  secondary  fuselage  vortex  and  the  stagnation  of  the  spanwise 
crossflow  component  on  the  wing.  Note  also  that  even  though  the  canard  vortex  had 
traveled  the  length  of  the  model,  the  pressure  gradient  around  it  was  approximately  equal  to 
that  around  the  wing  leading-edge  vortex.  This  further  indicates  that  it  was  the  energy  in  the 
canard  vortex  which  allowed  or  entrained  the  wing  vortex  to  reform.  The  -1.21  isobar 
marks  the  limit  of  the  probe  capabilities. 
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C.   CONCLUSION 

The  addition  of  the  canard  had  a  dramatic  effect  on  the  wing  flowfield  at  an  angle  of 
attack  of  22°.  Without  the  canard,  the  flow  across  the  wing  was  characterized  by  large 

areas  of  chaotic,  incoherent  flow  with  a  resultant  loss  in  lift.  The  addition  of  the  canard 
vortex  established  a  large  crossflow  component  moving  toward  the  wing  surface  which 
caused  the  flow  over  the  inboard  section  of  the  wing  to  reattach. 

In  addition  to  reattachment,  the  crossflow  provided  energy  to  the  wing  vortex  which 
delayed  its  breakdown.  Therefore,  with  the  canard  in  this  position  relative  to  the  wing,  the 
the  wing  leading-edge  vortex  was  enhanced  and  stabilized.  This  crossflow-induced 
breakdown  delay  substantiates  the  conclusions  made  by  Lacey,  in  Reference  2. 

The  energy  contained  in  the  canard  vortex  is  not  lost.  Consider  the  drag  force  as  a 
measure  of  the  energy  required  to  propel  two  models  with  the  same  wing  planform.  The 
first  model  employs  a  canard/wing  configuration  and  the  second  model  a 
wing/conventional-tail  configuration.  If  the  flow  over  the  wing  of  the  second  model 
remains  attached  (this  study  shows  it  would  not),  both  models  would  have  approximately 
the  same  drag.  The  difference  is  that  the  energy  (in  the  form  of  drag)  placed  in  the  tail 
vortex  of  the  second  model  would  be  lost  to  the  freestream.  In  the  canard/wing 
configuration,  some  of  the  energy  placed  in  the  canard  vortex  was  used  to  benefit  the  flow 
over  the  wing  thereby,  making  the  canard  configuration  more  efficient. 

Figure  31  is  provided  as  a  schematic  of  the  approximate  path  of  the  vortices  as  they 
moved  downstream  through  the  grids. 
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Freestream  Direction 

► 


Figure  31.  Vortex  Path 
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APPENDIX  A 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


DEF  SEG  :  CLEAR  ,  &HFE00 
GOTO  1000  '  User  program 
GOTO  900  '  Error  handling 
I  =  &HFE00'  Copyright  Hewlett-Packard 
PCIB.DIR$  =  ENVIRON$("PCIB") 
1$  =  PCIB.DIR$  +  "\PCIBILC.BLD" 
BLOAD  1$,  I 

CALL  I(PCIB.DIR$,  1%,  J%) 
IF  J%  =  0  THEN  GOTO  13 
PRINT  "Unable  to  load."; 


GOTO  4 'Begin  PCIB  Program  Shell 


1984, 1985 


PCIB.  SEG  =  I! 


PRINT  " 
END 

DEF  SEG  = 
I.C  =  20: 
GET. MEM  = 
PCIB.ERR$ 


(Error  #";  J! 


") 


10:  I.V  =  15 


17 

18  CALL  DEF. ERR (PCIB. ERR, 
PCIB.BASERR  =2  55 

19  ON  ERROR  GOTO  3 

20  J  =  -1 

21  1$  =  PCIB.DIR$  + 
2  2  CALL  O.S(I$) 

23  IF  PCIB. ERR  <>  0 

24  I  =  0 

2  5  CALL  I . V ( I , 
INITIALIZE . SYSTEM) 
26  IF  PCIB. ERR  <>  0 


PCIB. SEG:  O.S  =  5:  C.S  = 
L.P  =  25:  LD.FILE  =  30 

35:  L.S  =  40:  PANELS  =  45:  DEF. ERR  =  50 
=  STRING$(64,  32):  PCIB.NAME$  =  STRING$(16,  32) 
PCIB.ERR$,  PCIB.NAME$,  PCIB . GLBERR) : 


"\PCTB.SYN" 


THEN  ERROR  PCIB.BASERR 


READ. REGISTER, 


READ.SELFID, 


DEFINE, 


THEN  ERROR  PCIB.BASERR 


ENABLE. SYSTEM,   DISABLE . SYSTEM,   INITIALIZE, 


2  7   CALL  I.V  (I, 

POWER. ON) 

2  8  IF  PCIB. ERR  <>  0  THEN  ERROR  PCIB.BASERR 

2  9  CALL  I. V( I, MEASURE, OUTPUT, START, HALT) 

3  0  IF  PCIB. ERR  <>  0  THEN  ERROR  PCIB.BASERR 

31   CALL   I.V (I,    ENABLE. INT. TRIGGER,    DISABLE . INT . TRIGGER, 
ENABLE. OUTPUT,  DISABLE . OUTPUT) 
3  2  IF  PCIB. ERR  <>  0  THEN  ERROR 
3  3    CALL    I.V (I,    CHECK. DONE, 
SET. RANGE) 

<>  0  THEN  ERROR 

SET. MODE,  WRITE, 

<>  0  THEN  ERROR 

DELAY,  SAVE. SYSTEM,  J,  J) 

<>  0  THEN  ERROR  PCIB.BASERR 


34 
35 
36 
37 
38 
39 
40 
SET 


IF  PCIB. ERR 
CALL  I.V (I, 
IF  PCIB. ERR 
CALL  I.V (I, 
IF  PCIB. ERR 
1  =  1 

CALL    I • V ( I , 
SOURCE) 


PCIB.BASERR 
GET . STATUS , 

PCIB.BASERR 
CAL,  READ.CAL, 
PCIB.BASERR 


SET. FUNCTION 


STORE. CAL) 


SET.GATETIME 


SET. SAMPLES 


SET. SLOPE 
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41  IF  PCIB.ERR 
4  2  CALL  1.0(1, 

43  IF  PCIB.ERR 

44  CALL  I.C(I, 
4  5  IF  PCIB.ERR 
4  6  CALL  I.C(I, 

47  IF  PCIB.ERR 

48  CALL  1.0(1, 

49  IF  PCIB.ERR 

50  CALL  I.C(I, 

51  IF  PCIB.ERR 

52  I  =  2 

53  I  =  3 

54  CALL  I.V(I, 

55  IF  PCIB.ERR 

56  CALL  I.C(I, 

57  IF  PCIB.ERR 

58  CALL  I.C(I, 

59  IF  PCIB.ERR 

60  CALL  I.C(I, 

61  IF  PCIB.ERR 

62  CALL  I.C(I, 

63  IF  PCIB.ERR 

64  I  =  4 
65 


<>  0  THEN  ERROR 
FREQUENCY,  AUTO 
<>  0  THEN  ERROR 
INTERVAL,  RATIO 
<>  0  THEN  ERROR 
Rl,  RIO,  R100 


PCIB.BASERR 

FREQ,  PERIOD,  AUTO. PER) 

PCIB.BASERR 

TOTALIZE,  R100MILLI) 
PCIB.BASERR 
R1KIL0) 


<>  0  THEN  ERROR  PCIB.BASERR 
Rl OMEGA,  RIO OMEGA,  CHAN. A,  CHAN.B) 
<>  0  THEN  ERROR  PCIB.BASERR 
POSITIVE,  NEGATIVE,  COMN ,  SEPARATE) 
<>  0  THEN  ERROR  PCIB.BASERR 


ZERO. OHMS,  SET. SPEED,  J,  J) 

<>  0  THEN  ERROR  PCIB.BASERR 

DCVOLTS,  ACVOLTS,  OHMS,  R2  00MILLI) 

<>  0  THEN  ERROR  PCIB.BASERR 

R2,  R20,  R200,  R2KIL0) 

<>  0  THEN  ERROR  PCIB.BASERR 

R20KILO,  R200KILO,  R2MEGA,  R20MEGA) 

<>  0  THEN  ERROR  PCIB.BASERR 

AUTOM,  R2.5,  R12.5,  J) 

<>  0  THEN  ERROR  PCIB.BASERR 

SET. COMPLEMENT,   SET. DRIVER,   OUTPUT . NO . WAIT , 


CALL  I . V ( I , 
ENABLE . HANDSHAKE) 

66  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

67  CALL     I.V(I,      DISABLE. HANDSHAKE, 
SET. START. BIT,  SET . NUM. BITS) 

<>  0  THEN  ERROR  PCIB.BASERR 

SET. LOGIC 

<>  0  THEN 

POSITIVE, 

<>  0  THEN 

OC,  TTL, 

<>  0  THEN  ERROR  PCIB 

R2,  R3,  R4,  R5) 

<>  0  THEN  ERROR  PCIB 

R6,  R7,  R8,  R9) 

<>  0  THEN  ERROR  PCIB.BASERR 

RIO,  Rll, 

<>  0  THEN 

R14,  R15, 

<>  0  THEN 


SET. THRESHOLD, 


IF  PCIB.ERR 
CALL  I.V(I, 
IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
1  =  6 
CALL  I .  V  ( I 
SET. SYMMETRY) 
85  IF  PCIB.ERR 
8  6  CALL  I.V(I, 


68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 


SENSE,  J,  J,  J) 
ERROR  PCIB.BASERR 
NEGATIVE,  TWOS,  UNSIGNED) 
ERROR  PCIB.BASERR 
RO,  Rl) 

BASERR 

BASERR 


R12,  R13) 

ERROR  PCIB.BASERR 

R16,  J) 

ERROR  PCIB.BASERR 


SET. FREQUENCY,   SET . AMPLITUDE ,   SET. OFFSET, 


<>  0  THEN  ERROR  PCIB.BASERR 
SET. BURST. COUNT,  J,  J,  J) 


59 


87  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

88  CALL  I.C(I,  SINE,  SQUARE,  TRIANGLE,  CONTINUOUS) 

89  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

90  CALL  I.C(I,  GATED,  BURST,  J,  J) 

91  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
92 
93 


AUTOSCALE,    CALIBRATE,    SET . SENSITIVITY 


SET. TRIG. SLOPE 


GET. VERT. INFO 


CALC.WFTIME 


CALC. PERIOD, 


1  =  7 

CALL   I • V ( I , 
SET. VERT. OFFSET) 

94  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

95  CALL  I.V(I,   SET. COUPLING,   SET . POLARITY ,   SET . SWEEPSPEED, 
SET. DELAY) 

96  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

97  CALL     I.V(I,     SET. TRIG. SOURCE, 
SET. TRIG. LEVEL,  SET. TRIG. MODE) 

98  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

99  CALL   I.V(I,   GET. S INGLE. WF,   GET. TWO. WF, 
GET . TIMEBASE . INFO) 

100  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

101  CALL   I.V(I,   GET. TRIG. INFO,   CALC.WFVOLT, 
CALC.WF. STATS) 

102  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

103  CALL  I.V(I,   CALC.RISETIME,   CALC . FALLTIME , 
CALC. FREQUENCY) 

104  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

105  CALL     I.V(I,      CALC.PLUSWIDTH,      CALC . MINUSWIDTH 
CALC. OVERSHOOT,  CALC . PRESHOOT) 

106  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

107  CALL   I.V(I,   CALC.PK.TO.PK,   SET. TIMEOUT, 
MEASURE . S INGLE . WF ) 

108  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

109  CALL  I.V(I,  MEASURE . TWO . WF ,  J,  J,  J) 

110  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

111  CALL  1.0(1,  R10NANO,  R100NANO,  R1MICRO, 

112  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

113  CALL  1.0(1,  R100MICRO,  R1MILLI,  R10MILLI, 

114  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

115  CALL  I.C(I,  Rl,  RIO,  R20NANO,  R200NANO) 

116  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

117  CALL  I.C(I,  R2MICRO,  R20MICRO,  R200MICRO, 

118  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

119  CALL  I.C(I,  R20MILLI,  R200MILLI,  R2 ,  R20) 
12  0  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
121  CALL  I.C(I,  R50NANO,  R500NANO,  R5MICRO,  R50MICRO) 
12  2  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
123  CALL  I.C(I,  R500MICRO,  R5MILLI,  R50MILLI , 
12  4  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
12  5  CALL  I.C(I,  R5,  R50,  CHAN. A,  CHAN.B) 

IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

C(I,  EXTERNAL,  POSITIVE,  NEGATIVE,  AC) 


SCOPE. START 


126 

127  CALL  I 


R10MICRO) 
R100MILLI) 

R2MILLI) 


R500MILLI) 
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128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

800 

801 

899 

900 

901 

902 

903 

904 
ii  . 


IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
CALL  I.C(I, 
IF  PCIB.ERR 
1  =  8 

CALL  I.V(I, 
IF  PCIB.ERR 
CALL  C.S 
IF  PCIB.ERR 
1$  =  PCIB.DIR$ 
CALL  L.P(I$) 
IF  PCIB.ERR  <> 
1$  =  "DMM.01" : 


<>  0  THEN  ERROR  PCIB.BASERR 

DC,  TRIGGERED,  AUTO. TRIG,  AUTO. LEVEL) 

<>  0  THEN  ERROR  PCIB.BASERR 

XI,  XI 0,  STANDARD,  AVERAGE) 

<>  0  THEN  ERROR  PCIB.BASERR 


OPEN. CHANNEL,  CLOSE . CHANNEL,  J,  J) 
<>  0  THEN  ERROR  PCIB.BASERR 


<> 


0  THEN  ERROR  PCIB.BASERR 
+  "\PCIB.PLD" 


0  THEN  ERROR  PCIB.BASERR 

1  =  3:  J  =  0:  K  =  0:  L  = 
CALL  DEFINE (DMM.01,  1$,  I,  J,  K,  L) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1$  =  "Func.Gen.01" :  I  =  6:  J  =  0:  K  =  1: 
CALL  DEFINE (Func. Gen. 01,  1$,  I,  J,  K,  L) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1$  =  "Scope. 01" :  I=7:J=0:K=2:L 
CALL  DEFINE(Scope.01,  1$,  I,  J,  K,  L) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1$  =  "Counter. 01":  I=1:J=0:K=3: 
CALL  DEFINE (Counter. 01,  1$,  I,  J,  K,  L) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1$  =  "Dig. In. 01":  I=4:J=0:K=4:L=1 
CALL  DEFINE(Dig.In.01,  1$,  I,  J,  K,  L) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1$  =  "Dig. Out. 01":  I  =  4 :  J  =  1:  K  =  4:  L  =  1 
CALL  DEFINE (Dig. Out. 01,  1$, 
IF  PCIB.ERR  <>  0  THEN  ERROR 
1$  =  "Relay. Act. 01":  I=8:J=0:K=5:L=1 


L  =  1 


=  1 


L  =  1 


=  1:  K  =  4: 
I,  J,  K,  L) 
PCIB.BASERR 
J  =  0:  K  = 


CALL  DEFINE (RELAY. ACT. 01,  1$,  I,  J,  K,  L) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1$  =  "Relay. Mux. 01":  I  =  2:  J  =  0:  K  =  6: 
CALL  DEFINE (RELAY. MUX. 01,  1$,  I,  J,  K,  L) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1$  =  ENVIRON$( "PANELS")  +  "\PANELS.EXE" 
CALL  L.S(I$) 
GOTO  2 
IF  ERR  =  PCIB.BASERR  THEN  GOTO  903 


L  =  1 


"BASIC  error  #";  ERR; 


PRINT 
STOP 

TMPERR  =  PCIB.ERR:  IF  TMPERR  = 
PRINT  "PC  Instrument  error  #"; 
ERL 

905  PRINT  "Error:  " ;  PCIB.ERR$ 

906  IF   LEFT$ ( PCI B. NAME $,    1) 
"Instrument:  ";  PCIB.NAME$ 


occurred  in  line 


ERL 


0  THEN  TMPERR  =  PCIB.GLBERR 
TMPERR;  "  detected  at  line 


<>   CHR$(32)    THEN    PRINT 


61 


907  STOP 

908  COMMON  PCIB.DIR$,  PCIB.SEG 

9  09  COMMON  LD.FILE,  GET. MEM,  PANELS,  DEF.ERR 

910  COMMON  PCIB.BASERR,  PCIB.ERR,  PCIB.ERR$,  PCIB.NAME$, 
PCIB.GLBERR 

911  COMMON 
READ. REGISTER, READ . SELFID, DEFINE , INITIALIZE . SYSTEM, ENABLE . SY 
STEM , DISABLE . SYSTEM , INITIALIZE , POWER . ON , MEASURE , OUTPUT , START 
, HALT , ENABLE . INT . TRIGGER , DISABLE . INT . TRIGGER , ENABLE . OUTPUT , D 
I SABLE . OUTPUT , CHECK . DONE , GET . STATUS 

912  COMMON  SET. FUNCTION,  SET. RANGE,  SET. MODE,  WRITE . CAL 
READ.CAL,  STORE. CAL,  DELAY,  SAVE. SYSTEM,  SET.GATETIME 
SET. SAMPLES,  SET. SLOPE,  SET. SOURCE,  ZERO. OHMS,  SET. SPEED 
SET. COMPLEMENT,  SET. DRIVER,  OUTPUT. NO. WAIT,  ENABLE . HANDSHAKE 
DISABLE . HANDSHAKE 

913  COMMON  SET. THRESHOLD,  SET. START. BIT,  SET. NUM. BITS 
SET. LOGIC. SENSE,  SET. FREQUENCY,  SET. AMPLITUDE ,  SET. OFFSET 
SET. SYMMETRY,  SET . BURST. COUNT,  AUTOSCALE,  CALIBRATE 
SET. SENSITIVITY,  SET. VERT . OFFSET ,  SET . COUPLING ,  SET. POLARITY 
SET.SWEEPSPEED 

914  COMMON  SET. DELAY,  SET . TRIG. SOURCE  ,  SET . TRIG . SLOPE 
SET. TRIG. LEVEL,  SET . TRIG . MODE ,  GET . SINGLE . WF ,  GET . TWO . WF 
GET. VERT. INFO,  GET. TIMEBASE . INFO,  GET. TRIG . INFO ,  CALC.WFVOLT 
CALC.WFTIME,  CALC.WF. STATS ,  CALC. RISETIME,  CALC. FALLTIME 
CALC. PERIOD 

915  COMMON  CALC. FREQUENCY,  CALC . PLUSWIDTH  ,  CALC . MINUSWIDTH 
CALC. OVERSHOOT,  CALC . PRESHOOT,  CALC . PK. TO . PK,  SET. TIMEOUT 
SCOPE. START,  MEASURE . SINGLE . WF,  MEASURE . TWO . WF ,  OPEN. CHANNEL 
CLOSE. CHANNEL 

916  COMMON  FREQUENCY,  AUTO.FREQ,  PERIOD,  AUTO. PER,  INTERVAL 
RATIO,  TOTALIZE,  R100MILLI,  Rl ,  RIO,  R100,  R1KILO,  R10MEGA 
R100MEGA,  CHAN. A,  CHAN.B,  POSITIVE,  NEGATIVE,  COMN ,  SEPARATE 
DCVOLTS,  ACVOLTS,  OHMS,  R2  00MILLI,  R2 ,  R2  0,  R2  00,  R2KILO 
R20KILO  _ 

,  R200KILO 

917  COMMON  R2MEGA,  R20MEGA,  AUTOM,  R2 . 5 ,  R12.5,  POSITIVE 
NEGATIVE,  TWOS,  UNSIGNED,  OC ,  TTL,  RO ,  Rl,  R2  ,  R3  ,  R4  ,  R5 ,  R6 
R7,  R8,  R9,  RIO,  Rll,  R12,  R13 ,  R14 ,  R15,  R16,  SINE,  SQUARE 
TRIANGLE,  CONTINUOUS,  GATED,  BURST,  R10NANO,  R100NANO 
R1MICRO, 

R10MICRO,  R100MICRO 

918  COMMON  R1MILLI,  R10MILLI,  R100MILLI,  Rl ,  RIO,  R20NANO 
R200NANO,  R2MICRO,  R20MICRO,  R200MICRO,  R2MILLI,  R20MILLI 
R200MILLI,  R2,  R20,  R50NANO,  R500NANO,  R5MICRO,  R50MICRO 
R500MICRO,  R5MILLI,  R50MILLI,  R500MILLI,  R5 ,  R50,  CHAN . A 
CHAN . B , 

EXTERNAL,  POSITIVE 

919  COMMON  NEGATIVE,  AC,  DC,  TRIGGERED,  AUTO. TRIG,  AUTO. LEVEL, 
XI,  X10,  STANDARD,  AVERAGE 
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XPT(50) ,  YPT(50) ,  X(50) 


RELAY  MUX 


RELAY 


920   COMMON   DMM.01,    Func.Gen.01,   Scope. 01,    Counter. 01, 
Dig. In. 01,  Dig. Out. 01,  RELAY .ACT. 01 ,  RELAY. MUX. 01 

999  'End  PCIB  Program  Shell 

1000  REM  This  step  initialzes  the  HP  system 
1010  CLS 

1020  OPTION  BASE  1 

1030  DIM  P(5),  PA(50,  5),  PP(50,  5) 

Y(50),  YAW(50) 

1040  REM 

1050  CALL  INITIALIZE. SYSTEM (PGMSHEL. HPC) 

1060  REM 

1070  REM  SET  FUNCTIONON  THE   'DMM'   , 

ACTUATOR ' 

108  0  REM 

109  0  CALL  SET. FUNCTION (DMM.01,  DCVOLTS) 
1100  CALL  SET. RANGE (DMM.01,  AUTOM) 
1110  CALL  DISABLE. INT. TRIGGER (DMM.01) 

112  0  CALL  ENABLE. OUTPUT (RELAY. MUX. 01) 

113  0  CALL  ENABLE. OUTPUT (RELAY. ACT. 01) 
1140  REM  *************  PROGRAM  TRAVERSE 
1150  REM 

1160  REM        OPEN  THE  COM  PORT  AND  INITIALIZE 
SETTINGS 
1170  OPEN  "coml:1200,n,8,l,rs,cs,ds,cdn  FOR  RANDOM  AS  #1 

REM  SET  MOTOR  DEFAULT  VALUES 

DATA  2  000,2  000,2  000,2,2,2,0.00012  5, 0.00012  5, 0 

READ  VI,  V2,  V3,  Rl ,  R2 ,  R3 ,  CI,  C2 ,  C3 

REM  DEFINE  CHARACTERS  FOR  DATA  REDUCTION  ALGORITHM 

RN2$  =  "RENAME  AlRAW.DAT  " 

HEAD1$  =  "  #     X      Y 
P5       YAW   " 

FORMAT$  =  "##    ##.##  ##.## 
###.###   ###.###   ###.##" 
1250  PRINT 

********************************************* 

"**  USER  MUST  SELECT  'CAPS  LOCK7  FUNCTION  **" 
•I*******************************************" 

DISPLAY  MOTOR  DEFAULT  SETTINGS 


******************* 


THE  MOTOR 


1180 

1190 

1200 

1210 

1220 

1230 

P4 

1240 


PI 


P2 


000125 


P3 


###.###    ###.###    ###.### 


1260 
1270 
1280 
1290 
1300 
1310 


PRINT 

PRINT 

PRINT 

REM 

PRINT  " 

PRINT  " 
SETTINGS:11 
13  2  0  PRINT  " 
13  3  0  PRINT  " 
STEPS/SEC" 2) " 
1340  PRINT  " 
INCHES" 
1350  PRINT  " 
1360  PRINT 
1370  PRINT  "NOTE! ! 


************************************* 
INITIALIZED  VALUES  FOR  ALL  MOTOR 


VELOCITY  =  1000  STEPS/SEC" 
RAMP (MOTOR  ACCELERATION)  =  2 


(6000 


DEFAULT  INCREMENTAL  UNITS  ARE 


***********************************" 


USE  MANUAL  CONTROL  TO  INITIALIZE  PROBE 
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POSITION  BEFORE" 
1380  PRINT  " 


SELECTING  COMPUTER  CONTROLLED  MOVEMENT 


1390  PRINT 

14  00  INPUT  "MANUAL  CONTROL  OR  COMPUTER  CONTROL  (ENTER  'MAN7  or 

7CP7)" ;  CON$ 

1410  IF  CON$  =  "CP"  THEN  3490 

142  0  REM   OPTION  TO  CHANGE  DEFAULT  SETTINGS  OF  VELOCITY  OR 

ACCELERATION  RAMP 

14  3  0  PRINT 

14  4  0  PRINT 

14  50   PRINT   "   DO 

ACCELERATION  RAMP" 

1460  PRINT  " 

PRINT 

PRINT 


1470 

1480 

THE" 

1490  PRINT 

1500  PRINT 

FOR  ANY  MOTOR . " 

1510  PRINT 

1520  PRINT 

1530  PRINT 

1540  INPUT 

(Y  or  N) "; 

1550  IF  D$ 

1560  IF  D$ 

1570  REM 

1580  REM 

1590  PRINT 

1600  PRINT 

1610  INPUT 


"IF  7NO' 


YOU   WANT   TO   CHANGE   THE   VELOCITY   OR 


DEFAULT  SETTINGS?  (Y  or  N) " 


THIS  PROGRAM  WILL  THEN  LET  YOU  DEFINE 


"DISTANCE  YOU  WANT  TO  MOVE  (IN  INCHES).  IF  'YES7," 
"YOU  CAN  CHANGE  ANY  OR  ALL  OF  THE  DEFAULT  SETTINGS 


"DO  YOU  WANT  TO  CHANGE  ANY  OF  THE  DEFAULT  SETTINGS? 

D$ 

_  hY"  THEN  1590 

=  "N"  THEN  2220 

****  OPERATOR  SELECTED  MOTOR  VARIABLES  ***** 


"WHICH  DEFAULT 
FOR  ACCEL  RAMP)";  L 
1620  ON  L  GOTO  1690,  1930 
163  0  PRINT  "DO  YOU  WANT  TO 
N)  " 

1640  INPUT 
1650  IF  V$ 
1660  PRINT 
RAMP?  (Y  or 
1670 
1680 
1690 
1700 
1710 
or  3) 
1720 
1730 
1740 


VALUE?  (ENTER  ' 1 ' FOR  VELOC  OR  72 


CHANGE  THE  DEFAULT  VELOCITY?  (Y  OR 


V$ 

=  "Y"  THEN  1690 
"DO  YOU  WANT  TO 
N)  " 

THEN  1990 
THEN  1450 


CHANGE  THE  DEFAULT  ACCELERATION 


=  "v" 
=  nNn 


IF  R$ 
IF  R$ 
PRINT 
PRINT 
INPUT  "WHICH  MOTOR  VELOCITY  DO  YOU  WISH  TO  CHANGE?  (1,2, 

";  J 

ON  J  GOTO 

PRINT 

PRINT 


1730,  1830,  1880 


64 


1750  INPUT 

1760  PRINT 

1770  PRINT 

1780  PRINT 

1790  PRINT 

(Y  OR  N) " 

1800  INPUT  V$ 

1810  IF  V$  =  "Y"  THEN 

1820  IF  V$  =  "N"  THEN 

183  0  PRINT 

1840  PRINT 

18  50  INPUT 

1860  PRINT 

1870  GOTO  1780 

1880  PRINT 

1890  PRINT 

1900  INPUT  "ENTER 

1910  PRINT 

1920  GOTO  1780 

193  0  PRINT 

1940  PRINT 

1950  INPUT  "WHICH 

2 ,  or  3 ) " ;  K 


"ENTER  DESIRED  VELOCITY  OF  MOTOR  #1";  VI 


"DO  YOU  WANT  TO  CHANGE  VELOCITY  OF  ANOTHER  MOTOR? 


1690 
1430 


"ENTER  DESIRED  VELOCITY  OF  MOTOR  2";  V2 


DESIRED  VELOCITY  OF  MOTOR  #3";  V3 


MOTOR  ACCEL  RAMP  DO  YOU  WANT  TO  CHANGE?  (1, 


"ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  #1";  Rl 


"DO  YOU  WANT  TO  CHANGE  THE  ACCEL  RAMP  OF  ANOTHER 


THEN 
THEN 


1930 
1450 


ON  K  GOTO  1970,  2060,  2120 

PRINT 

PRINT 

INPUT 

PRINT 

PRINT 

PRINT 
MOTOR?  (Y  or  N) ?" 
2  03  0  INPUT  RM$ 

IF  RM$  =  "Y" 

IF  RM$  =  "N" 

PRINT 

PRINT 

INPUT 

PRINT 

PRINT 

GOTO  2000 

PRINT 

PRINT 

INPUT 

PRINT 

PRINT 

GOTO 

REM 

REM  DEFINE  DISTANCE  TO  MOVE  MOTOR 

PRINT 


1960 
1970 
1980 
1990 
2000 
2010 
2020 


2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 


"ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  #2";  R2 


"ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  #3";  R3 


2000 
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2  210  PRINT 
2220  PRINT 
2230  REM  INITIALIZE  MOTOR  INCREMENTS  TO  ZERO 


2240 
2250 
2260 


11  = 

12  = 

13  = 


2270  PRINT 


8 


0 


N 


DEFINE  WHICH  MOTOR  YOU  WANT  TO  MOVE 

NOTE!!!   A  POSITIVE  ('+')  INCREMENT  TO 
MOVES  TRAVERSER  AWAY  FROM  THAT  PARTICULAR 


********************************************************* 

2290  PRINT  "  ** 
**» 

2300  PRINT  "  ** 
**» 

2  310  PRINT  "  ** 
A  MOTOR    **" 
2320  PRINT  "  ** 
MOTOR     **" 
2330  PRINT  "  ** 
**»» 

2  340  PRINT  "  **   —  MOTOR  #1  MOVES  THE  PROBE  UPSTREAM  AGAINST 

THE  FLOW   **" 

2350  PRINT  "  **   —  MOTOR  #2  MOVES  THE  PROBE  TOWARD  THE  ACCESS 

WINDOW    **" 

2360  PRINT  "   **    —  MOTOR  #3  MOVES  THE  PROBE  VERTICALLY 

DOWNWARD  **" 

2       3       7       0PRINT" 

************************************************************** 

2  380  PRINT 

2  39  0  PRINT 

2400  INPUT  "WHICH  MOTOR  DO  YOU  WANT  TO  MOVE?  (1,2,  or  3 ) " ;  L 

2410  ON  L  GOTO  2420,  2680,  2970 

242  0  PRINT 

24  3  0  PRINT 

2  44  0  PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #1?" 

2450  PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********" 

2460  INPUT  II 

2470  PRINT 

2480  PRINT  "  *********************************" 

2490  PRINT 

2500  PRINT  "SUMMARY  OF  OPERATOR  INPUTS:" 

2510  PRINT  "         MOTOR  #1    VELOCITY  =  " ;  VI 

2520  PRINT  "  ACCELERATION  RAMP  =  ";  Rl 

2530  PRINT  "  INCREMENTAL  DISTANCE  =  ";  II; 

"INCHES" 

2540  PRINT  "  *********************************" 

2  550  PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or 

N)  " 

2560  PRINT 

2  57  0  PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.   ENTER  'Y'  TO 
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RETURN" 
2580  PRINT 


2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

"INCHES" 

2  810  PRINT 


"TO 

V$ 

=  Hyil 

3410 


THEN  14  3  0 


"DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N) 

C$ 

—  ii  y» 

=  »N" 


THEN 
THEN 


2220 
3260 


VARIABLE  SELECTION  SUBROUTINE 
INPUT 
IF  V$ 
GOSUB 
PRINT 
PRINT 
INPUT 
IF  C$ 
IF  C$ 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 

REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 
PRINT  "  *********************************" 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


911 


"HOW  FAR  DO 
"  ********* 

12 


YOU  WANT  TO  MOVE  MOTOR  #2?" 

(ENTER  DISTANCE  IN  INCHES)  ********** 


"SUMMARY 
ii 

ti 

ii 


OF  OPERATOR  INPUTS:" 

MOTOR  #2    VELOCITY  =  ";  V2 

ACCELERATION  RAMP  =  " ;  R2 
INCREMENTAL  DISTANCE  =  "; 


12 


*********************************  ii 


2820 

2830 

2840 

N)  " 

2850 

2860 


PRINT 
PRINT 
PRINT 

PRINT 
PRINT 


RETURN" 
287  0  PRINT 


2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 


INPUT 
IF  V$ 
GOSUB 
PRINT 
PRINT 
INPUT 
IF  C$ 
IF  C$ 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 


"DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or 


"ENTER  'N'  TO  START  MOTOR  MOVEMENT.   ENTER  'Y'  TO 

"TO  VARIABLE  SELECTION  SUBROUTINE." 

V$ 

=  "Y"  THEN  1430 

3410 

"DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N)?" 

C$ 

=  "Y"  THEN  2220 

=  "N"  THEN  3260 


"HOW  FAR  DO 
••  ********* 

13 


YOU  WANT  TO  MOVE  MOTOR  #3?" 

(ENTER  DISTANCE  IN  INCHES)  *********" 
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REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 
PRINT  "  *****************************#***» 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


3030 

3040 

3050 

3060 

3070 

3080 

3090 

"INCHES" 

3100  PRINT 


"SUMMARY  OF  OPERATOR  INPUTS:" 


MOTOR  #3 


VELOCITY  =  ";  V3 
ACCELERATION  RAMP  =  ' 
INCREMENTAL  DISTANCE    " ;  13; 


R3 
ii  . 


3110 

3120 

3130 

3140 

N)" 

3150 

3160 


PRINT 
PRINT 
PRINT 
PRINT 

PRINT 
PRINT 


*********************************** 


"DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or 


"ENTER  'N'  TO  START  MOTOR  MOVEMENT.   ENTER 


TO 


RETURN" 
3170  PRINT 


"TO 

V$ 

=  ii  Y" 

3410 


THEN  14  3  0 


VARIABLE  SELECTION  SUBROUTINE 
INPUT 
IF  V$ 
GOSUB 
PRINT 
PRINT 
INPUT  "DO  YOU  WANT  TO  INPUT  ANOTHER  MANUAL  MOTOR  MOVEMENT 

M$ 

=  "Y"  THEN  2210 


"DO  YOU  WANT  TO  INPUT  COMPUTER  CONTROLLED  MOTOR 


3180 

3190 

3200 

3210 

3220 

3230 

(Y  or  N) "; 

3240  IF  M$ 

3250  PRINT 

3  2  60  PRINT 

MOVEMENT?" 

3  27  0  PRINT 

3  28  0  PRINT 

ZEROIZED." 

3290  PRINT  "PROGAM  WILL  LET  YOU  CHOOSE  MANUAL  OR  CP-CONTROLLED 

MOVEMENT . " 

3  3  00  PRINT 

3  310  PRINT 

3  32  0  INPUT 

or  N) ";  N$ 

3330  IF  N$ 


*********  NOTE!!!  *********  ii 
ALL  PREVIOUS  MOTOR  INCREMENT  INPUTS 


HAVE  BEEN 


******  (if  'NO',  THE  PROGRAM  WILL  END).  *****" 
"DO  YOU  WANT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT  (Y 
=  "Y"  THEN  3500 


3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
END 
REM 
PRINT 


****************************" 

THE  PROGRAM  HAS  ENDED." 
****************************** 


*******  MOTOR  MOVEMENT  SUBROUTINE  ********* 
#1,  "&":  PRINT  #1,  "E";  "Cl=" ;  CI;  ":C2  =  " 


C2; 


":C3=";  C3 
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3430  PRINT  #1,  "11=";  II;  ":V1=";  VI;  ":R1=";  Rl; 
3440  PRINT  #1,  ":I2=" ;  12;  ":V2=";  V2 ;  " :R2=";  R2 
3450  PRINT  #1,  "13=";  13;  ":V3=";  V3 ;  ":R3=";  R3 ;  ":@" 

3  4  60  RETURN 

347  0  REM    ******************************************** 

3  480  REM    ******************************************** 

3  4  90  PRINT 

3  500  REM  *******  COMPUTER  CONTROLLED  MOVEMENT  ******* 

3  510  PRINT 

3  52  0  PRINT  "THE  PRESSURE  DATA  WILL  BE  WRITTEN  TO  FILES  ON 

DRIVE  'A'  " 

3  53  0  PRINT 

3540  PRINT  "YOU  WILL  BE  ASKED  TO  INPUT  FILE  NAMES  FOR  THESE." 

3  550  PRINT 

3560  INPUT  "IS  A  FORMATTED  DISK  IN  DRIVE  'A'?   PRESS  'ENTER7 

TO  CONTINUE";  D$ 

3570  PRINT 

3580  PRINT 

3  59  0  PRINT 


***************************************** 
**  NOTE  ! ! !  **' 

**  COMPUTER  CONTROLLED  MOVEMENT  **' 
**        IS  PROGRAMMED  WITH  A  **' 

**  DEFAULTED  NEGATIVE  MOTOR  INCREMENT  **' 
**  (i.e.  MOTOR  #3  WILL  MOVE  UPWARD  **' 
**  BY  ENTERING  A  (+)  DISTANCE) .  **' 
***************************************** 


3600  PRINT 

3610  PRINT 

3  62  0  PRINT 

3  63  0  PRINT 

3  64  0  PRINT 

3  650  PRINT 

3  660  PRINT 

3  67  0  PRINT 

3680  PRINT 

3  690  REM  SET  INITIAL  MOVEMENT  DISTANCE  AND  NUMBER  OF  DATA 

POINTS  TO  ZERO 

3700  HT  =  0 

3710  WD  =  0 

3720  DIST  =  0 

3730  XPT  =  0 

3740  YPT  =  0 

3750  N  =  0 

3760  PRINT 

3770  PRINT 

3780  INPUT  "WHAT  IS  THE  DIMENSION  (  X  ,  Y  )   (IN  INCHES)  THAT 

YOU  WANT  TO  MEASURE.";  WD,  HT 

3  79  0  PRINT 

3  8  00  INPUT  "WHAT  IS  THE  STEP  (IN  INCHES)  THAT  YOU  WANT  TO 

MOVE.";  DIST 

3810  YPT  =  INT(HT  /  DIST)  +  1 

3820  XPT  =  INT (WD  /  DIST)  +  1 

383  0  N  =  XPT  *  YPT 

3840  PRINT 

3850  PRINT  "THERE  ARE   ";  XPT;  "  *  ";  YPT;  "  =  " ;  N;  "  POINTS 

TO  BE  MEASURED  " 
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3860 

3870 

3880 

3890 

3900 

3910 

3920 

NAMES 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 


PRINT 

INPUT  "ARE  THE  NUMBER 

IF  C$  =  "N"  THEN  3780 

CLS 

N  =  XPT 

IF  (N  <  1)  OR  (N  > 

REM  ***  GENERATING 


OF  POINTS  IS  OK. (Y  OR  N) 


C$ 


99)  GOTO  3780 
STRING  STRING 


SEGMENTS  FOR  DATA  FILE 


MID$(STR$(1) 
MID$(STR$(N) 
"XXXXXX" 
DAT" 


** 
** 


STRING 
ENDING 


NUMBER 
STRING 


it  -^ii 
NUMBER  "N" 


"DATA  FILES  WILL  BE  INCREMENTED  FROM:" 


(X$  +  E$  +  EX$) 


"ENTER 
;  F2$ 


DATA  FILE  NAME  (6  CHARACTERS  MAX  —  NO 


>  6  OR  LEN(F2$)  <  1  GOTO  4  03  0 


B$  =  MID$(STR$(1) ,  2):  REM 

E$  =  MID$(STR$(N) ,  2):  REM 

X$  = 

EX$  = 

CLS 

PRINT 

PRINT 

PRINT  (X$  +  B$  +  EX$) ;  "   To 

PRINT 

PRINT 

INPUT 
EXTENSION) " 
4  04  0  PRINT 
4050  PRINT 
4060  IF  LEN(F2$) 
4070  CLS 

4080  PRINT  N;  "DATA  FILES  WILL  BE  GENERATED  AND  INCREMENTED  AS 
FOLLOWS : " 
4  09  0  PRINT 

PRINT 

PRINT  (F2$  +  B$  +  EX$) ;  "   To   " ;  (F2$  +  E$  +  EX$) 

PRINT 

PRINT 

INPUT  "ARE  THE  NUMBER  OF  POINTS  AND  FILE  NAMES  OK. (Y  OR 

C$ 

IF  C$  =  "N"  GOTO  3780 

IF  C$  =  "Y"  GOTO  4180 

GOTO  414  0 

CLS 

PRINT 

PRINT 

REM  SET  INITIAL  POSITION  DATA 

X(l)  =  -DIST 

Y(l)  =  -DIST 

FOR  IX  =  2  TO 

X(IX)  =  0 

NEXT  IX 

FOR  JY  =  2 

Y(JY)  =  0 

NEXT  JY 

FOR  I  =  1  TO  XPT 


4100 
4110 
4120 
4130 
4140 
N)  "  ; 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 


XPT  +  1 


TO  YPT  +  1 
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4302  II  =  0 

4304  12  =  0 

4306  13  =  0 

4  310  FOR  J  =  1  TO  YPT 

4  320  REM   MOTOR   CP-CONTROLLED  MOTOR  MOVEMENT 

4330  II  =  0 

4340  12  =  0 

4350  13  =  0 

4  3  60  REM  EACH  POINT  TAKE  10  TIMES  READINGS 

4370  X(I  +  1)  =  X(I)  +  DIST 

4380  XPT(J)  =  X(I  +  1) 

4390  Y(J  +  1)  =  Y(J)  +  DIST 

4400  YPT(J)  =  Y(J  +  1) 

4405  INPUT  "  ADJUST  THE  WHEEL  TO  MAKE  THE  P2  =P3 ,  INPUT  THE  YAW 

ANGLE";  YAW (J) 

4408  PRINT 

4410  INPUT  "  PRESS  'ENTER'  TO  START  THE  MEASUREMENT";  MOVE$ 

4420  REM 

443  0  REM  READ  FIVE  CHANNELS  AND  DISPLAY  THE  DATA 

4440  REM 

4450  STEPPER  =  4 

4460  SWITCH  =  3 

4470  HOMER  =  8 

4480  DELAY1  =  .1 

4490  DELAY2  =  1 

4500  REM  SET  THE  S.V  PORT  TO  #4 

4510  FOR  IL  =  1  TO  3 

4520  THYME  =  TIMER 

4530  CALL  OUTPUT (RELAY . ACT . 01 , STEPPER) 

4  54  0  CHKTIME  =  TIMER 

4550  IF  CHKTIME  <  (THYME  +  DELAY1)  GOTO  4540 

4560  CALL  OPEN . CHANNEL (RELAY .ACT . 01 ,  SWITCH) 

4570  CLS 

4580  NEXT  IL 

4  590  PRINT 

4600  PRINT  "  NOW  IS  POINT  ";  J 

4  610  REM  START  MEASURE  FROM  PORT  4  TO  PORT  8 

4  62  0  FOR  JJ  =  1  TO  5 

4  63  0  CALL  OUTPUT (RELAY . ACT . 01 , STEPPER) 

4  64  0  CHKTIME  =  TIMER 

4  650  IF  CHKTIME  <  (THYME  +  DELAY2)  GOTO  4  64  0 

4  660  REM  EACH  PORT  SAMPLE  10  TIMES 

4670  FOR  II  =  1  TO  10 

4  68  0  ROUT  =  1 

4  69  0  CALL  OUTPUT (RELAY . MUX . 01 , ROUT) 

4700  CALL  MEASURE (DMM. 01,  VOLTS) 

4  710  PA (II,  JJ)  =  VOLTS 

4720  NEXT  II 

4730  CALL  OPEN . CHANNEL (RELAY . ACT . 01 ,  SWITCH) 
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4740  IF  JJ  =  5  THEN  4760 

4750  NEXT  JJ 

4760  REM  HOME  THE  S.V  PORT  TO  #48 

4  77  0  CALL  OUTPUT (RELAY . ACT . 01 , HOMER) 

478  0  CALL  OPEN . CHANNEL (RELAY . ACT . 01 ,  HOMER) 

4790  REM 

48  00  REM  DISPLAY  THE  SAMPLE  DATA 

4810  REM 

4820  PRINT  HEAD1$ 

4830  FOR  IS=  1  TO  10 

4840     PRINT     USING 

FORMAT$;IS,XPT(J) ,YPT(J) , PA (IS , 1) , PA (IS , 2 ) , PA (IS , 3 ) , PA ( IS , 4 ) 

, PA (IS ,5) ,YAW(J) 

4850  NEXT  IS 

4860  REM 

487  0  REM  AVERAGE  THE  DATA 

4880  REM 

4890  FOR  JA  =  1  TO  5 

4  9  00  TOTAL  =  0 

4910  FOR  IA  =  1  TO  10 

4920  TOTAL  =  TOTAL  +  PA(IA,  JA) 

4930  NEXT  IA 

494  0  AVERAGE  =  TOTAL  /  10 

4950  P(JA)  =  AVERAGE 

4960  NEXT  JA 

4970  PRINT 

498  0  PRINT  "THE  AVERAGES  ARE:  " 

5000  PRINT  HEAD1$ 

5010  FOR  JD  =  1  TO  5 

5020  PP(J,  JD)  =  P(JD) 

5030  NEXT  JD 

5040  PRINT  USING  FORMAT $ ;  J;  XPT(J);  YPT(J);  PP(J,  1);  PP(J, 

2);  PP(J,  3);  PP(J,  4);  PP(J,  5);  YAW(J) 

5045  PRINT 

5048  PRINT  USING  "THE  NULLING  ERROR  IS  +#.####";  PP(J,  3)  - 
PP(J,  2) 

5049  PRINT 

5050  PRINT  "DO  YOU  WANT  RE-MEASURE  AGAIN  (Y  /  N) " 
5060  PRINT 

5062  PRINT  "IF  'Y'  WILL  RE-SAMPLE  AGAIN." 

5064  PRINT 

5070  INPUT  "IF  'N'  WILL  MOVE  THE  TRAVERSER  STEP  UPWARD  (WAIT 

7  SEC  ) " ;  C$ 

5075  PRINT 

5080  IF  C$  =  "Y"  THEN  4405 

5082  IF  C$  =  "N"  THEN  5090 

5084  GOTO  5070 

5090  IF  J  =  YPT  THEN  5160 

5100  REM 
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5110  REM  MOVE  THE  TRAVERSER  STEP  UPWARD. 

5120  REM 

5130  13  =  -DIST 

5140  GOSUB  3410 

5150  NEXT  J 

5160  REM***  STORE  DATA  BEFORE  NEXT  SAMPLE*** 

5170  OPEN  "A:\RAW.DAT"  FOR  OUTPUT  AS  #2 

5180  PRINT  #2,  HEAD1$ 

5190  FOR  ID  =  1  TO  YPT 

5200  PRINT  #2,  USING  FORMAT $ ;  ID;  XPT(ID);  YPT(ID);  PP(ID,  1)  ; 

PP(ID,  2);  PP(ID,  3);  PP(ID,  4);  PP(ID,  5);  YAW(ID) 

5210  NEXT  ID 

5220  CLOSE  #2 

5230  REM  ***  GENERATING  INCREMENTED  DATA  FILE  NAME 

5240  IF  (I  >  10)  OR  (I  =  10)  THEN  1$  =  MID$ (STR$ ( I ) ,  2) 

5250  IF  (I  <  10)  THEN  1$  =  (MID$ (STR$ (0) ,  2)  +  MID$ (STR$ (I) , 

2)) 

5260  FI2$  =  (F2$  +  1$  +  EX$) 

5270  PRINT 

5280  PRINT  "  WRITING  DATA  FILE  " ;  FI2$ 

5290  DF2$  =  RN2$  +  FI2$ 

53  00  REM  **  RENAME   DATA  FILE 

5310  SHELL  DF2$ 

5320  REM 

53  3  0  REM  MOVE  THE  TRAVERSER  TO  THE  NEXT  SAMPLE  POSITION 

5340  REM 

5350  PRINT 

5360  IF  I  =  XPT  THEN  5430 

53  70  INPUT  "THEN  PRESS  'ENTER'  FOR  NEXT  COLUMN  SAMPLE (  9  0  SEC) 

11 ;  MOVE$ 

5390  12  =  -DIST 

5400  13  =  HT 

5410  GOSUB  3410 

5420  NEXT  I 

5430  CLS 

544  0  PRINT  "ALL  MOVEMENTS  COMPLETE" 

5450  PRINT 

5460  PRINT 

547  0  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER 

MOVEMENT  (Y  OR  N)?" 

5480  PRINT 

5490  PRINT  "IF   'Y',   THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL 

CONTROL  SUBROUTINE." 

5500  PRINT  "IF  'N' ,  THE  PROGRAM  WILL  END." 

5510  PRINT 

5520  INPUT  "ANOTHER  MOVEMENT";  R$ 

5530  IF  R$  =  "Y"  THEN  1370 

5540  IF  R$  =  "N"  THEN  3370 
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DEF  SEG  :  CLEAR  ,  &HFE00 

GOTO  1000  '  User  program 

GOTO  900  '  Error  handling 

I  =  &HFE00'  Copyright  Hewlett-Packard 

PCIB.DIR$  =  ENVIRON$("PCIB") 

1$  =  PCIB.DIR$  +  "\PCIBILC.BLD" 

BLOAD  1$,  I 

CALL  I(PCIB.DIR$,  1%,  J%):  PCIB.SEG  = 

IF  J%  =  0  THEN  GOTO  13 

PRINT  "Unable  to  load."; 

PRINT  " 

END 


APPENDIX  B 


GOTO  4 'Begin  PCIB  Program  Shell 


1984, 1985 


1% 


(Error  #";  J%;  ") 


DEF  SEG  = 
I.C  =  20: 
GET. MEM  = 
PCIB.ERR$ 


10:  I.V  =  15 


PCIB.SEG:  O.S  =  5:  C.S  = 

L.P  =  25:  LD.FILE  =  30 

35:  L.S  =  40:  PANELS  =  45:  DEF. ERR  =  50 

=  STRING$(64,  32):  PCIB.NAME$  =  STRING$(16, 


32) 


17 

18  CALL  DEF. ERR (PCIB. ERR,  PCIB.ERR$,  PCIB.NAME$,  PCIB.GLBERR) 
PCIB.BASERR  =2  55 

19  ON  ERROR  GOTO  3 

20  J  =  -1 

21  1$  =  PCIB.DIR$  +  "\PCIB.SYN" 
2  2  CALL  O.S(I$) 

2  3  IF  PCIB. ERR  <> 
24  I  =  0 

2  5  CALL  I . V ( I , 
INITIALIZE . SYSTEM) 
26  IF  PCIB. ERR  <> 


0  THEN  ERROR  PCIB.BASERR 


READ. REGISTER,    READ.SELFID, 


DEFINE, 


0  THEN  ERROR  PCIB.BASERR 


27   CALL  I.V(I,   ENABLE. SYSTEM,   DISABLE . SYSTEM,   INITIALIZE, 

POWER. ON) 

2  8  IF  PCIB. ERR  <>  0  THEN  ERROR  PCIB.BASERR 

2  9  CALL  I.V (I, MEASURE, OUTPUT, START, HALT) 

3  0  IF  PCIB. ERR  <>  0  THEN  ERROR  PCIB.BASERR 

31  CALL   I.V(I,    ENABLE. INT. TRIGGER,    DISABLE 
ENABLE. OUTPUT,  DISABLE . OUTPUT) 

32  IF  PCIB. ERR  <>  0  THEN  ERROR 
3  3  CALL  I.V (I,  CHECK. DONE 
SET. RANGE) 

<>  0  THEN  ERROR 

SET. MODE,  WRITE 

<>  0  THEN  ERROR 

DELAY,  SAVE. SYSTEM,  J,  J) 

<>  0  THEN  ERROR  PCIB.BASERR 


34 
35 
36 
37 
38 
39 
40 


IF  PCIB. ERR 
CALL  I.V(I, 
IF  PCIB. ERR 
CALL  I.V(I, 
IF  PCIB. ERR 
1  =  1 

CALL    I .  V  ( I 


PCIB.BASERR 
GET. STATUS, 

PCIB.BASERR 
CAL,  READ.CAL, 
PCIB.BASERR 


INT. TRIGGER, 


SET. FUNCTION 


STORE. CAL) 


SET.GATETIME 


SET. SAMPLES 


SET. SLOPE 


SET. SOURCE) 
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43 

44 


41  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
4  2  CALL  I.C(I,  FREQUENCY,  AUTO . FREQ ,  PERIOD,  AUTO. PER) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
CALL  I.C(I,  INTERVAL,  RATIO,  TOTALIZE,  R100MILLI) 
4  5  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

46  CALL  I.C(I,  Rl,  RIO,  R100,  R1KILO) 

47  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

48  CALL  I.C(I,  R10MEGA,  R100MEGA,  CHAN. A,  CHAN.B) 

49  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

50  CALL  I.C(I,  POSITIVE,  NEGATIVE,  COMN,  SEPARATE) 

51  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 


52 
53 


I  = 
I  = 


63 
64 
65 


SET. DRIVER,   OUTPUT . NO . WAIT, 


SET. THRESHOLD, 


54  CALL  I.V(I,  ZERO. OHMS,  SET. SPEED,  J,  J) 

55  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

56  CALL  1.0(1,  DCVOLTS,  ACVOLTS ,  OHMS,  R2  00MILLI) 

57  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

58  CALL  I.C(I,  R2,  R20,  R200,  R2KIL0) 

59  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

60  CALL  I.C(I,  R20KILO,  R200KILO,  R2MEGA,  R20MEGA) 

61  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

62  CALL  I.C(I,  AUTOM,  R2 . 5 ,  R12.5,  J) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1  =  4 

CALL  I.V(I,   SET. COMPLEMENT, 
ENABLE . HANDSHAKE ) 

66  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

67  CALL     I.V(I,      DISABLE. HANDSHAKE, 
SET. START. BIT,  SET. NUM. BITS) 

68  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

69  CALL  I.V(I,  SET. LOGIC. SENSE,  J,  J,  J) 

70  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

71  CALL  1.0(1,  POSITIVE,  NEGATIVE,  TWOS,  UNSIGNED) 

72  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

73  CALL  I.C(I,  OC,  TTL,  RO ,  Rl) 

74  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

7  5  CALL  I.C(I,  R2,  R3 ,  R4 ,  R5) 

76"  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

77  CALL  I.C(I,  R6,  R7 ,  R8 ,  R9) 

78  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

79  CALL  I.C(I,  RIO,  Rll,  R12,  R13) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

I,  R14,  R15,  R16,  J) 
IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 
1  =  6 
CALL   I.V(I,   SET. FREQUENCY,   SET . AMPLITUDE , 
SET. SYMMETRY) 
85  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

8  6  CALL  I.V(I,  SET. BURST. COUNT,  J,  J,  J) 


80 

81  CALL  I 

82 

83 

84 


SET. OFFSET, 
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87  IF  PCI B. ERR  <>  0  THEN  ERROR  PCIB.BASERR 

88  CALL  1.0(1,  SINE,  SQUARE,  TRIANGLE,  CONTINUOUS) 

89  IF  PCI B. ERR  <>  0  THEN  ERROR  PCIB.BASERR 

90  CALL  I.C(I,  GATED,  BURST,  J,  J) 

91  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

92  I  =  7 

9  3    CALL   I.V(I,    AUTOSCALE,    CALIBRATE,    SET . SENSITIVITY , 
SET. VERT. OFFSET) 

94  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

95  CALL  I.V(I,   SET. COUPLING,   SET . POLARITY ,   SET . SWEEPSPEED , 
SET. DELAY) 

96  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

97  CALL     I.V(I,     SET. TRIG. SOURCE,     SET . TRIG . SLOPE , 
SET. TRIG. LEVEL,  SET. TRIG. MODE) 

98  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

99  CALL   I.V(I,   GET. SINGLE. WF,   GET. TWO. WF,   GET . VERT . INFO , 
GET . TIMEBASE . INFO) 

100  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

101  CALL   I.V(I,   GET. TRIG. INFO,   CALC.WFVOLT,   CALC.WFTIME, 
CALC.WF. STATS) 

102  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

103  CALL  I.V(I,   CALC.RISETIME,   CALC . FALLTIME ,   CALC. PERIOD, 
CALC. FREQUENCY) 

104  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

105  CALL     I-V(I,     CALC.PLUSWIDTH,      CALC . MINUSWIDTH , 
CALC. OVERSHOOT,  CALC . PRESHOOT) 

106  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

107  CALL   I.V(I,   CALC.PK.TO.PK,   SET. TIMEOUT,   SCOPE. START, 
MEASURE . S INGLE . WF ) 

108  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

109  CALL  I.V(I,  MEASURE . TWO . WF ,  J,  J,  J) 

110  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

111  CALL  I.C(I,  R10NANO,  R100NANO,  R1MICRO,  R10MICRO) 

112  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

113  CALL  I.C(I,  R100MICRO,  R1MILLI,  R10MILLI,  R100MILLI) 

114  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

115  CALL  I.C(I,  Rl,  RIO,  R2  0NANO,  R2  00NANO) 

116  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

117  CALL  I.C(I,  R2MICRO,  R20MICRO,  R200MICRO,  R2MILLI) 

118  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

119  CALL  I.C(I,  R2  0MILLI,  R2  00MILLI,  R2 ,  R2  0) 
12  0  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

121  CALL  I.C(I,  R50NANO,  R500NANO,  R5MICRO,  R50MICRO) 

12  2  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

12  3  CALL  I.C(I,  R500MICRO,  R5MILLI,  R50MILLI,  R500MILLI) 

12  4  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

12  5  CALL  I.C(I,  R5,  R50,  CHAN. A,  CHAN.B) 

126  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

127  CALL  I.C(I,  EXTERNAL,  POSITIVE,  NEGATIVE,  AC) 
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128  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

129  CALL  I.C(I,  DC,  TRIGGERED,  AUTO. TRIG,  AUTO. LEVEL) 
13  0  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

131  CALL  I.C(I,  XI,  X10,  STANDARD,  AVERAGE) 
13  2  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

133  1=8 

134  CALL  I.V(I,  OPEN. CHANNEL,  CLOSE . CHANNEL,  J,  J) 

135  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

13  6  CALL  C.S 

137  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

138  1$  =  PCIB.DIR$  +  "\PCIB.PLD" 

139  CALL  L.P(I$) 

140  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

141  1$  =  "DMM.01":  I=3:J=0:K=0:L=1 

142  CALL  DEFINE (DMM.01,  1$,  I,  J,  K,  L) 

14  3  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

144  1$  =  "Func.Gen.01":  I  =  6:J  =  0:K=1:L=1 

145  CALL  DEFINE (FUNC. GEN. 01,  1$,  I,  J,  K,  L) 
14  6  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

147  1$  =  "Scope. 01":  1=7:  J  =  0:  K=2:  L=l 

148  CALL  DEFINE(SCOPE.01,  1$,  I,  J,  K,  L) 

149  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

150  1$  =  "Counter. 01":  I=1:J=0:K=3:L=1 

151  CALL  DEFINE (COUNTER. 01,  1$,  I,  J,  K,  L) 

152  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

153  1$  =  "Dig. In. 01":  I=4:J=0:K=4:L=1 

154  CALL  DEFINE(DIG.IN.01,  1$,  I,  J,  K,  L) 

155  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

156  1$  =  "Dig. Out. 01":  I=4:J=1:K=4:L=1 

157  CALL  DEFINE(DIG.OUT.01,  1$,  I,  J,  K,  L) 

158  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

159  1$  =  "Relay. Act. 01":  I=8:J=0:K=5:L=1 

160  CALL  DEFINE (RELAY. ACT. 01,  1$,  I,  J,  K,  L) 

161  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

162  1$  =  "Relay. Mux. 01":  I=2:J=0:K=6:L=1 

163  CALL  DEFINE (RELAY. MUX. 01,  1$,  I,  J,  K,  L) 

164  IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB.BASERR 

800  1$  =  ENVIRON$( "PANELS")  +  "\PANELS.EXE" 

801  CALL  L.S(I$) 
899  GOTO  2 

9  00  IF  ERR  =  PCIB.BASERR  THEN  GOTO  903 

901  PRINT  "BASIC  error  #";  ERR;  "  occurred  in  line  " ;  ERL 

902  STOP 

903  TMPERR  =  PCIB.ERR:  IF  TMPERR  =  0  THEN  TMPERR  =  PCIB.GLBERR 

904  PRINT  "PC  Instrument  error  #";  TMPERR;  "  detected  at  line 
"  ;  ERL 

905  PRINT  "Error:  ";  PCIB.ERR$ 

906  IF   LEFT$  ( PCI  B.  NAME  $,    1)    <>   CHR$(32)    THEN   PRINT 
"Instrument:  ";  PCIB.NAME$ 
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907  STOP 

908  COMMON  PCIB.DIR$,  PCIB.SEG 

909  COMMON  LD.FILE,  GET. MEM,  PANELS,  DEF.ERR 

910  COMMON  PCIB.BASERR,  PCIB.ERR,  PCIB.ERR$,  PCIB.NAME$, 
PCIB.GLBERR 

911  COMMON 
READ . REGISTER , READ . SELFID , DEFINE , INITIALI ZE . SYSTEM , ENABLE . S Y 
STEM , DISABLE . SYSTEM , INITIALI ZE , POWER . ON , MEASURE , OUTPUT , START 
, HALT , ENABLE . INT . TRIGGER , DISABLE . INT . TRIGGER , ENABLE . OUTPUT , D 
I SABLE . OUTPUT , CHECK . DONE , GET . STATUS 

912  COMMON  SET. FUNCTION,  SET. RANGE,  SET. MODE,  WRITE. CAL 
READ.CAL,  STORE. CAL,  DELAY,  SAVE. SYSTEM,  SET.GATETIME 
SET. SAMPLES,  SET. SLOPE,  SET. SOURCE,  ZERO. OHMS,  SET. SPEED 
SET. COMPLEMENT,  SET. DRIVER,  OUTPUT . NO .WAIT,  ENABLE . HANDSHAKE 
DISABLE . HANDSHAKE 

913  COMMON  SET. THRESHOLD,  SET . START . BIT,  SET. NUM. BITS 
SET. LOGIC. SENSE,  SET . FREQUENCY ,  SET .AMPLITUDE ,  SET. OFFSET 
SET. SYMMETRY,  SET. BURST. COUNT,  AUTOSCALE,  CALIBRATE 
SET. SENSITIVITY,  SET . VERT . OFFSET ,  SET . COUPLING,  SET. POLARITY 
SET.SWEEPSPEED 

914  COMMON  SET. DELAY,  SET. TRIG . SOURCE ,  SET . TRIG . SLOPE 
SET. TRIG. LEVEL,  SET. TRIG. MODE ,  GET. SINGLE . WF,  GET. TWO. WF 
GET. VERT. INFO,  GET. TIMEBASE. INFO,  GET. TRIG. INFO,  CALC.WFVOLT 
CALC.WFTIME,  CALC . WF . STATS ,  CALC . RISETIME ,  CALC . FALLTIME 
CALC. PERIOD 

915  COMMON  CALC . FREQUENCY ,  CALC . PLUSWIDTH ,  CALC . MINUSWIDTH 
CALC. OVERSHOOT,  CALC. PRESHOOT,  CALC. PK. TO. PK,  SET. TIMEOUT 
SCOPE. START,  MEASURE . SINGLE . WF ,  MEASURE . TWO . WF ,  OPEN. CHANNEL 
CLOSE. CHANNEL 

916  COMMON  FREQUENCY,  AUTO.FREQ,  PERIOD,  AUTO. PER,  INTERVAL 
RATIO,  TOTALIZE,  R100MILLI,  Rl ,  RIO,  R100,  R1KILO,  R10MEGA 
R100MEGA,  CHAN. A,  CHAN.B,  POSITIVE,  NEGATIVE,  COMN,  SEPARATE 
DCVOLTS,  ACVOLTS,  OHMS,  R2  00MILLI,  R2 ,  R20,  R2  00,  R2KILO 
R2  0KILO  _ 

,  R200KILO 

917  COMMON  R2MEGA,  R2 OMEGA,  AUTOM,  R2 . 5 ,  R12.5,  POSITIVE 
NEGATIVE,  TWOS,  UNSIGNED,  OC,  TTL,  RO ,  Rl ,  R2 ,  R3 ,  R4 ,  R5 ,  R6 
R7,  R8,  R9,  RIO,  Rll,  R12 ,  R13,  R14 ,  R15,  R16,  SINE,  SQUARE 
TRIANGLE,  CONTINUOUS,  GATED,  BURST,  R10NANO,  R100NANO 
R1MICRO,   _ 

R10MICRO,  R100MICRO 

918  COMMON  R1MILLI,  R10MILLI,  R100MILLI,  Rl ,  RIO,  R20NANO 
R2  00NANO,  R2MICRO,  R2  0MICRO,  R2  00MICRO,  R2MILLI ,  R2  0MILLI 
R200MILLI,  R2,  R20,  R50NANO,  R500NANO,  R5MICRO,  R50MICRO 
R500MICRO,  R5MILLI,  R50MILLI,  R500MILLI,  R5 ,  R50,  CHAN . A 
CHAN . B ,   _ 

EXTERNAL,  POSITIVE 

919  COMMON  NEGATIVE,  AC,  DC,  TRIGGERED,  AUTO. TRIG,  AUTO. LEVEL, 
XI,  X10,  STANDARD,  AVERAGE 
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92  0   COMMON   DMM.01,    FUNC.GEN.01,    SCOPE. 01,    COUNTER. 01, 
DIG. IN. 01,  DIG. OUT. 01,  RELAY . ACT. 01 ,  RELAY. MUX. 01 

999  'End  PCIB  Program  Shell 

1000  REM 

1010  REM  This  step  initialzes  the  HP  system 

1020  CLS 

1030  OPTION  BASE  1 

1040  DIM  P(10),  PA(50,  6),  PP(50,  6),  XPT(40),  CAL(40) 

1050  CALL  INITIALIZE. SYSTEM ( PGMSHEL. HPC) 

1060  REM 

1070  REM  All  PC  devices  now  have  an  initial  state 

1080  REM  Set  function  on  the  DMM  and  Relay  MUX 

1090  REM 

1100  CALL  SET. FUNCTION (DMM. 01,  DCVOLTS) 

1110  CALL  SET. RANGE (DMM. 01,  AUTOM) 

112  0  CALL  DISABLE. INT. TRIGGER (DMM. 01) 

113  0  CALL  ENABLE. OUTPUT (RELAY. MUX. 01) 

1140   FORMAT$  =   "##    ##.####    ##.####    ##.####    ##.#### 

1200  FOR  I  =  1  TO  10 

1210  CAL(I)  =  0! 

12  2  0  NEXT  I 

1510  REM 

1520  REM  READ  THE  VOLTAGE  OF  48TH   CHANNEL  AND  DISPLAY  THE 

DATA 

1530  REM 

154  0  PRINT  "  CHOOSE  6  POINTS" 

1550  PRINT 

1550  PRINT  "THE  CALIBRATION  WILL  BE  STORES  IN  'CAL.DAT'" 

1560  REM 

1570  REM  Begin  sampling  loop 

1580  REM 

1600  FOR  J  =  1  TO  1 

1610  PRINT 

1630  FOR  JJ  =  1  TO  6 

1631  INPUT  "INPUT  THE  CALIBRATION  PRESSURE";  CAL(JJ) 

1632  INPUT  "PRESS  'ENTER'  TO  START   MEASUREMENT";  MOVE$ 
1640  FOR  II  =  1  TO  10 

1650  ROUT  =  1 

1660  CALL  OUTPUT (RE LAY. MUX. 01, ROUT) 

167  0  CALL  MEASURE (DMM. 01,  VOLTS) 

1680  PA(II,  JJ)  =  VOLTS 

1690  NEXT  II 

1700  IF  JJ  =  6  THEN  1740 

1730  NEXT  JJ 

1740  REM 

1750  REM  DISPLAY  THE  SAMPLE  DATA 

1760  REM 

1780  FOR  IS=  1  TO  10 
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1790     PRINT     USING 

FORMAT$;IS,PA(IS,l) , PA (IS , 2 ) , PA (IS , 3 ) , PA ( IS , 4 ) , PA ( IS , 5) , PA (I 

S,6) 

1800  NEXT  IS 

1810  REM 

18  2  0  REM  AVERAGE  THE  DATA 

18  3  0  REM 

184  0  FOR  JA  =  1  TO  6 

18  50  TOTAL  =  0 

1860  FOR  IA  =  1  TO  10 

187  0  TOTAL  =  TOTAL  +  PA(IA,  JA) 

1880  NEXT  IA 

189  0  AVERAGE  =  TOTAL  /  10 

1900  P(JA)  =  AVERAGE 

1920  NEXT  JA 

19  3  0  PRINT 

194  0  PRINT  "THE  AVERAGE  ARE:  " 

2000  FOR  JD  =  1  TO  6 

2010  PP(J,  JD)  =  P(JD) 

2020  NEXT  JD 

2055  PRINT  USING  FORMAT $ ;  J;  PP(J,  1);  PP(J,  2);  PP(J,  3); 

PP(J,  4) ;  PP(J,  5) ;  PP(J,  6) 

2  07  0  PRINT 

2080  INPUT  "DO  YOU  WANT  RE-MEASURE  AGAIN  ?  (Y  /  N)";  C$ 

2090  IF  C$  =  "Y"  THEN  1580 

2101  REM***  STORE  DATA  BEFORE  NEXT  SAMPLE*** 

2102  OPEN  "A: \CAL.DAT"  FOR  OUTPUT  AS  #2 

2106  FOR  ID  =  1  TO  6 

2107  PRINT  #2,  USING  FORMAT $ ;  ID;  PP(J,  ID);  CAL(ID) 

2108  NEXT  ID 

2109  CLOSE  #2 
2210  NEXT  J 
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APPENDIX  C 


$storage: 2 

$debug 

************************************************************ 

***** 

*  THIS   PROGRAM   CONVERTS   THE   VOLTAGE   OF   TRANSDUCER   INTO 
PHYSICAL  * 

*  PRESSURE, VELOCITY, YAW  ANGLE  AND  PITCH  ANGLE.  THOSE  DATA  ARE 
* 

*  USED   FOR  PLOT  PROGRAM  LATER. 
* 

************************************************************ 
***** 

CHARACTER* 14  FNAME 

CHARACTER* 14  NAME 

CHARACTER* 14  OUTFILE 

CHARACTER* 2  A (50) 

CHARACTER* 80  ST 

REAL  K,INTR 

INTEGER  COLS , RWS , DTPTS 

DATA  A/'01' , '02' , '03' , '04' , '05' ,  /06 

*  '10' , 'll' , '12' , 

*  '  19 ' 

*  '28' 

*  '37' 

*  '46' 
WRITE    (*, ' (A\) ') 

(LEFT/RIGHT)  =  ' 
READ 
WRITE 


'29'  ,  '30' ,  '31' 
'38'  ,  '39' ,  '40' 
'47'  ,  '48' ,  '49' 


'14', '15' 


'07' , '08' , '09' 
'16'  '17'  '18' 


# 


'32'  ,  '33' 
'41' , '42' 
'50'/ 
OF   COLUMNS 


'34'  ,  '35'  ,  '36' 
'43'  ,  '44'  ,  '45' 


IN   THE   GRID 


(UP/ DOWN)  =  ' 
READ 
WRITE 
READ 
WRITE 

READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
curve  [cm/v] 


*,*)  COLS 
(*/'(a\)')   '   #   OF   DATA   POINTS 


IN  A   COLUMN 


* , * )  RWS 

(*,'(A\)')  '  DATA  FILE  NAME?  (A: FILEXXXX. EXT) :  ' 

*, ' (A14) ')  NAME 

(*,'(A\)')  '  Initial  ambient  pressure  [in.  Hg] : 


*,*)  PI 
(*,'(A\)') 
*,*)  PF 
(*,'(A\)') 
*  *)  TI 

(*,'(A\)') 
*,*)  TF 

(*,'(A\)')  '  Tunnel  calibration  factor,  K: 
*,  *)  K 
(*,'(A\)') 


Final  ambient  pressure  [in.  Hg 
Initial  temperature  [deg  F] :  ' 
Final  temperature  [deg  F] :  ' 


Slope  from  pressure  calibration 
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READ  (*,*)  SLOPE 

WRITE  (*, '  (A\) ')  '  Intercept  from  pressure  calibration: 

READ  (*,*)  INTR 

WRITE  (*,'(A\)')  '  Tunnel  delta-p  for  test-section  q: 

READ  (*,*)  QM1FAC 

*  COMPUTE  THE  AVERAGE  ATMOSPHERIC  PRESSURE  GIVEN  inHg  CONVERT 
TO  psf 

PATM=(PI+PF) *35.3631 

R=1716.5 
q    ******  "E"  is  blockage  correction  ************** 

E=0.0123 
c  ******  COMPUTE  TEST  SECTION  AVERAGE  TEMPERATURE 

T=(TI+TF)/2.+460 
C  ******  TEST  SECTION  DENSITY 

RO=PATM/(R*T) 

DTPTS=RWS*COLS 

*  OPEN  A  NEW  FILE  TO  STORE  THE  REDUCED  DATA 

write (*,' (A\) ') '  Output  data  file?  [A: FILENAME. EXT] :  ' 
read(*, ' (al4) ')  outfile 
OPEN ( 2 , FILE=outf ile , STATUS= ' new ' ) 
WRITE  (2,900)  DTPTS 
900    FORMAT(I5) 

*  OPEN  A  SEQUENTIAL  OF  DATA  FILE  *  BE  SURE  THE  FILE  IN  PPROBE 
HAS  6  ELEMENTS 

DO  2  0  1=1, COLS 

NAME (9: 10) =A (I) 
FNAME=NAME 
OPEN ( 1 , FILE=FNAME) 
READ(1, 100,END=2  0)ST 
100     FORMAT(A65) 

15       READ(1, 1000, END=30)NO,X,Y, VI, V2,V3,V4,V5, BETA 
1000     FORMAT ( 12 , F7 . 2 , F6 . 2 , 5F9 . 3 , F8 . 2 ) 

*  CONVERT   THE   VOLTAGE   TO   PRESSURE   IN   LBF/FT**2       ** 
lcm/H2O=2.04  61  PSF  *** 

Pl=calvp (VI , SLOPE , INTR) *2 . 0461+PATM 
P2=calvp(V2, SLOPE, INTR) *2. 0461+PATM 
P3=calvp(V3, SLOPE, INTR) *2. 0461+PATM 
P4=calvp (V4 , SLOPE , INTR) *2 . 0461+PATM 
P5=calvp(V5, SLOPE, INTR) *2 . 0461+PATM 

IF  ( (Pl-(P2+P3)/2)  .LT.  0.0)  THEN 

*  AN   ATTEMP   TO   ACCOUNT   FOR   THE   DELTA   P   NOT   EQUAL   TO 
Z  ERO  ************** 

*  P1=ABS(P1) 

*  P3=ABS(P3) 

*  P2=ABS(P2) 
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*  CALCULATE  THE  PITCH  ANGLE  IN  DEGREES 

P=(P4-P5)/(P1-P2) 
ALPHA=FPITCH(P) 

*  CALCULATE  THE  VELOCITY  IN  FT/SEC 

YSLOP=FYSLOP (ALPHA) 

VELM=SQRT( (2*YSLOP* (ABS ( (P1-P3) ) ) )/ (RO*K) ) 

VEL=VELM*(1+E) 

*  CALCULATE  THE  LOCAL  DYNAMIC  PRESSURE 

QM1=QM1FAC*2 . 0461/K 
QM=RO*VEL**2/2. 
Q1=QM1*(1+2*E) 
Q=QM*(1+2*E) 

*  CALCULATE  THE  YAW  ANGLE  IN  DEGREES 

c  **********  betaO  is  tunnel  cal  correction  ***** 
betaO  =   4 . 

YAW=FYAW ( BETA  -  betaO) 
c 
c  *****   Calculate  the  velocity  components 

betar  =  yaw*0. 017453 
C  *****   alphaO  is  tunnel  cal  correction  ***** 

alphaO  =  -2. 

pitch   =  alpha-alphaO 

alphar  =  pitch*0 . 017453 

vely  =  vel*sin(alphar) 

velx  =  vel*cos (alphar) *sin(betar) 
c 

*  CALCULATE  THE  TOTAL  PRESSURE  IN  LBF/IN**2 

PTC=FPT (ALPHA) 

PT1=P1-Q*PTC 

PT=PT1/14  4. 

CPT= ( PT1-PATM-Q1) /Ql 

IF  (CPT  .LT.  -3.0)  THEN 

CPT=-3.0 
ENDIF 

*  CALCULATE  THE  STATIC  PRESSURE  IN  LBF/IN**2 

PS1=PT1-Q 
PS=PS1/144. 
CPS= ( PS1-PATM) /Ql 
GO  TO  1190 
ELSE 

*  CALCULATE  THE  PITCH  ANGLE  IN  DEGREES 

P=(P4-P5)/(P1-P2) 
ALPHA=FPITCH(P) 

*  CALCULATE  THE  VELOCITY  IN  FT/SEC 

YSLOP=FYSLOP (ALPHA) 

VELM=SQRT( (2*YSLOP* ( ( (P1-P2) ) ) )/ (RO*K) ) 
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VEL=VELM*(1+E) 

*  CALCULATE  THE  LOCAL  DYNAMIC  PRESSURE 

QM1=QM1FAC*2 . 0461/K 
QM=RO*VEL**2/2. 
Q1=QM1*(1+2*E) 
Q=QM*(1+2*E) 

*  CALCULATE  THE  YAW  ANGLE  IN  DEGREES 

c  **********  betaO  is  tunnel  cal  correction  ***** 
betaO  =   4. 

YAW=F YAW ( BETA  -  betaO) 
c 
c  *****   Calculate  the  velocity  components 

betar  =  yaw*0. 017453 
C  *****   alphaO  is  tunnel  cal  correction  ***** 

alphaO  =  -2. 

pitch  =  alpha-alphaO 

alphar  =  pitch*0. 017453 

vely  =  vel*sin(alphar) 

velx  =  vel*cos (alphar) *sin (betar) 
c 

*  CALCULATE  THE  TOTAL  PRESSURE  IN  LBF/IN**2 

PTC=FPT (ALPHA) 

PT1=P1-Q*PTC 

PT=PT1/14  4. 

CPT= ( PT1-PATM-Q1 ) /Ql 

*  CALCULATE  THE  STATIC  PRESSURE  IN  LBF/IN**2 

PS1=PT1-Q 
PS=PS1/144. 
CPS=(PS1-PATM)/Q1 
ENDIF 

119  0 

WRITE (2, 2  000) X , Y , VEL, velx, vely , YAW, pitch , pt , cpt , ps , cps 
2000     FORMAT (11F10. 4) 

GO  TO  15 
3  0     CLOSE (1) 
2  0     CONTINUE 
CLOSE (2) 
STOP 
END 
************************************************************* 

*  THIS  FUNCTION  CONVERTS  THE  VOLTAGE  TO  PHYSICAL  PRESSURE 

FUNCTION  calvp(X, SLOPE, INTR) 
REAL  INTR 

calvp=X*SLOPE+INTR 
END 
************************************************************* 

*  THIS  FUNCTION   CALCULATES  THE  PITCH  ANGLE 

FUNCTION  FPITCH(X) 

FPITCH=3.7  59+5  3.7  568*X-1.3  08  5*X**2-1.658  3*X**3 
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*  -0.8061*X**4+16.5115*X**5 

END 
************************************************************ 

*  THIS  FUNCTION  CALCULATES  THE  VELOCITY  PRESSURE  COEFFICENT 

FUNCTION  FYSLOP(X) 
IF(X.LT.-10)THEN 

FYSLOP=0.981-0.0102*X-3.000E-4*X**2-2.500E-6*X**3 
ELSE  IF( (X.GE.-10) .AND. (X.LE.10) ) THEN 

FYSLOP=0.98-0.006*X+2.000E-4*X**2 
ELSE 

FYSLOP=0.9801-0.003  5*X-1.14  3E-4*X**2+5.8  3  3E-6*X**3 
END  IF 
END 
************************************************************* 

*  THIS  FUNCTION  CALCULATES  THE  YAW  ANGLE 

FUNCTION  FYAW(X) 
c********  NOTE:   180  deg  indicates  zero  yaw 
FYAW=180-x 

END 
************************************************************* 

*  THIS  FUNCTION  CALCULATES  THE  TOTAL  PRESSURE  COEFFICIENT 

FUNCTION  FPT(X) 
IF(X.LE.-30)  THEN 

FPT=-0.01 
ELSE  IF( (X.GT.-30) .AND. (X.LT.-20) )  THEN 

FPT=0.02+1.00E-3*X 
ELSE  IF( (X.GE.-20) .AND. (X.LE.30) )  THEN 

FPT=0 
ELSE 

FPT=0.03-1.00E-3*X 
END  IF 
END 
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